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A ductile-to-brittle  transition  can  occur  in  polymers  as  a result  of  many  factors 
due  to  their  unique  characteristics.  Physical  aging,  or  structural  relaxation,  is  a time-and- 
temperature-dependent,  as  well  as,  thermoreversible  phenomenon  that  occurs  in 
amorphous  materials  below  the  glass  transition,  i.e.,  T . Physical  aging  can  induce  a 
ductile-to-brittle  transition  in  amorphous  thermoplastics. 

The  objective  of  this  dissertation  is  to  study  the  ductile-to-brittle  transition  in 
amorphous  thermoplastics  due  to  molar  mass  and  thermal  history  by  using  a quantitative 
ffactography  technique  at  a macroscopic  and  a molecular  level.  The  influence  of  end 
groups  on  controlling  the  physical  aging  effect  is  also  studied.  Polycarbonate  with  two 
different  end  groups  and  different  ranges  of  molar  mass  were  used  in  the  study. 
Polycarbonates  were  characterized  using  gel  permeation  chromatography  (GPC), 
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differential  scanning  calorimetry  (DSC),  dynamic  mechanical  spectroscopy  (DMS),  and 
density  measurement.  Quantitative  fractography  at  the  macroscopic  level  was  done  by 
using  optical  microscopy  to  measure  the  sizes  of  different  regions  on  fracture  surfaces  of 
unnotched  samples  which  were  caused  to  fail  under  tension.  Fracture  toughness  (K1C) 
was  calculated  from  the  relationship  between  these  sizes  and  fracture  stress.  Quantitative 
fractography  at  the  molecular  level  was  done  by  atomic  force  microscopy  (AFM).  Fractal 
dimension  was  calculated  from  the  AFM  images.  Physical  aging  was  conducted  at  T - 
40°C  for  seven  days  for  the  DMS  experiments  and  at  Tg-20°C  for  seven  days  for  density 
measurements  and  tensile  tests. 

The  size  of  end  groups  effected  the  Tg,  the  shape  of  the  secondary  relaxation,  the 
activation  energy  of  the  secondary  relaxation  and  primary  relaxation  at  the  Tg,  and  the 
density  of  polycarbonate.  Large  end  groups  increase  Tg  and  decrease  the  density  of 
polycarbonate  at  the  same  range  of  molar  mass.  Large  end  groups  induce  broader 
secondary  relaxation  with  lower  activation  energy  but  higher  activation  energy  at  Tg.  The 
rates  of  aging  determined  from  the  increase  of  storage  modulus  as  a function  of  aging 
time  for  both  end  groups  are  not  significantly  different.  A ductile-to-brittle  transition  as  a 
function  of  molar  mass  can  be  observed  from  tensile  properties,  fracture  toughness,  and 
fractography.  A ductile-to-brittle  transition  as  a function  of  thermal  history  can  be 
observed  only  in  high  molar  mass  polycarbonate  from  tensile  properties  and  fracture 
toughness. 
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CHAPTER  1 
INTRODUCTION 


Amorphous  thermoplastic  polymers  have  been  used  extensively  because  of  their 
versatility.  Their  applications  range  from  household  goods  to  more  technically  advanced 
products  used  in  the  electronics  and  aerospace  industries.  Extensive  research  has  been 
undertaken  in  order  to  better  utilize  their  unique  properties  and/or  to  design  the  molecular 
structures  corresponding  to  their  applications. 

Fracture  behavior  is  one  of  the  important  characteristics  to  be  understood  prior  to 
the  design  process.  Polymers  display  very  complicated  fracture  behavior  due  to  their 
viscoelastic  nature.  Polymers  with  the  same  structure  can  exhibit  ductile  or  brittle 
fractures,  depending  on  their  thermal  history.  Understanding  the  effect  of  internal 
variables,  e.g.,  chemical  structures,  and  external  variables,  e.g.,  thermal  history,  on 
fracture  behavior  can  assist  a great  deal  in  the  design  process. 

Fractography  is  one  of  the  techniques  applied  to  study  fracture  behavior.  It 
involves  studying  the  detailed  morphology  of  fracture  surfaces  of  materials  at  the 
macroscopic,  microscopic,  or  molecular  level.  It  is  possible  to  identify  the  origin  of  a 
fracture,  its  direction  and  speed  of  propagation,  the  fracture  stress,  and  the  mode  of 
fracture  from  the  appearance  of  the  fracture  surface.  Therefore,  fractography  is  very 
useful  in  post-failure  analysis  to  solve  problems  in  the  field  of  material  engineering. 
Fractography  of  polymers  at  the  macroscopic  and  the  microscopic  level  has  been  studied 
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extensively  using  different  tools  such  as  the  optical  microscope  and  the  scanning  electron 
microscope.  However,  little  has  been  done  at  the  molecular  level.  Limited  publications 
are  available  that  attempt  to  quantify  the  information  obtained  from  fractography  or  to 
correlate  it  to  polymer  properties. 

The  invention  of  the  scanning  tunneling  microscope  and  the  atomic  force 
microscope  has  enabled  the  study  of  the  morphology  of  materials  at  the  nanometer  level. 
The  application  of  these  tools  is  widespread  in  several  areas  of  material  research,  e.g.,  in 
semiconductors,  surface  modification,  and  biomaterials.  In  the  field  of  fracture  research, 
these  tools  can  elucidate  the  fracture  mechanism  at  the  molecular  level,  which  may  be 
useful  in  molecular  design.  Furthermore,  by  using  these  tools  to  get  the  images  of 
fracture  surfaces  at  the  molecular  level,  fractal  dimension,  a quantitative  parameter  to 
describe  fracture  patterns,  can  be  determined. 

The  objective  of  this  dissertation  is  to  study  the  fracture  behavior  of  amorphous 
thermoplastic  polymers.  Comparisons  will  be  made  between  different  internal  variables, 
i.e.,  different  end  groups  and  different  molar  masses,  and  between  different  external 
variables,  i.e.,  thermal  histories.  Ductile-to-brittle  transition  caused  by  both  internal  and 
external  variables  will  be  investigated.  Fractography  at  the  macroscopic  level  using  the 
optical  microscope  and  at  the  molecular  level  using  the  atomic  force  microscope  will  be 
studied.  Fracture  toughness  (KIC)  and  fractal  dimension  will  be  determined  using 
quantitative  fractography  at  the  macroscopic  level  and  molecular  level,  respectively.  The 
relationship  between  fracture  toughness  and  fractal  dimension  will  be  studied  as  well. 


CHAPTER  2 
BACKGROUND 


Fracture  Behavior  of  Polymers 

According  to  the  theoretical  strength-of-solid  concept,  which  assumes  that  a 
material  is  a perfectly  elastic  isotropic  solid,  the  maximum  tensile  stress  (atheo)  will  be 
about  one-tenth  of  Young’s  modulus  (E)  (Kinloch  and  Young,  1983).  In  reality,  the 
maximum  strength  of  materials  is  usually  lower  due  to  either  inherent  or  processing 
flaws.  They  are  generally  of  the  order  of  E/50  to  E/100  for  isotropic  polymers  as  shown 
in  Table  2.1.  Therefore,  it  will  be  very  important  to  understand  the  fracture  behavior  of 
materials  and  the  relationship  between  strength  and  flaws. 

The  fracture  behavior  of  polymers  is  of  great  interest  in  the  field  of  engineering. 
Numerous  attempts  to  understand  the  mechanisms  and  to  predict  the  behavior,  as  well  as 
to  analyze  the  causes  of  failure,  have  been  made.  Polymers  exhibits  complicated  fracture 
behavior.  There  are  many  factors  that  influence  the  fracture  behavior  of  polymers,  e.g., 
molar  mass,  viscoelasticity,  and  type  of  polymers. 

Effect  of  Molar  Mass 

Molar  mass  is  a crucial  parameter  for  polymers.  To  possess  good  mechanical 
properties,  polymers  need  the  molar  mass  to  at  least  equal  the  Mc  (critical  molar  mass)  so 
the  polymers  will  have  sufficient  entanglement  density.  The  relationship  between  log 
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Table  2.1  Comparison  of  the  values  of  Young’s  modulus,  measured  tensile  strength,  and 
theoretical  tensile  strength  of  different  polymers 


Material 

E (MPa) 

CTtheo  (MPa) 

Of  (MPa) 

Typical  glassy  polymer 
(polymethylmethacrylate) 

3000 

300 

50 

Typical  semicrystalline  polymer 
(high  density  polyethylene) 

-2000 

-200 

20 

Typical  thermosetting  polymer 
(epoxy  resin) 

3500 

350 

70 

Typical  polymer  fiber 
(nylon  6) 

60000 

6000 

2000 

Source:  Kinloch  and  Young  (1983) 


melt  viscosity  and  log  molar  mass  is  universal.  The  slope  of  the  plot  changes  from  1 to 
3.4  at  Mc.  Fracture  behavior  also  depends  on  molar  mass;  generally,  high  molar  mass 
polymers  are  more  ductile  than  low  molar  mass  polymers  assuming  other  variables  are 
the  same. 

Mudrich  and  Wilchester  (1994)  studied  the  effect  of  molar  mass  on  the  fracture 
toughness  and  the  crack  propagation  rate  of  polycarbonate  (PC).  PC  with  Mw  of  18,  21, 
23,  and  26  kg/mole  were  used.  Double  torsion  geometry  was  used  in  the  study  in  order  to 
create  brittle  fracture  and  to  apply  the  linear  elastic  fracture  mechanics  (LEFM)  concept. 
The  results  showed  that  fracture  toughness  increases  with  molar  mass.  A marked 
increase  was  found  after  the  molar  mass  of  21  kg/mole  or  above  a critical  molar  mass. 
Crack  propagation  rates  were  found  to  be  comparable  at  the  molar  masses  of  21,  23,  and 
26  kg/mole.  They  are  all  lower  than  the  rate  at  the  molar  mass  of  1 8 kg/mole.  Again,  a 
marked  decrease  in  crack  propagation  rate  was  found  above  a critical  molar  mass. 

Pitman  and  Ward  (1979)  determined  fracture  toughness  of  different  PC  with  a 
wide  range  of  molar  masses  by  using  compact  tension  geometry.  The  results  show  that 
fracture  toughness  increases  as  a function  of  molar  mass.  There  was  a dramatic  increase 
in  fracture  toughness  between  Mn  5000  and  8000  g/mole. 
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Hui  and  Kramer  (1995)  proposed  the  micromechanics  model,  which  is  based  on 
the  cross-tie  fibrils  in  the  craze  microstructures.  These  cross-tie  fibrils  act  as  lateral  load 
bearers  and  they  transfer  stress  between  main  fibrils.  By  using  this  model,  the  fracture 
toughness  (Gic)  is  shown  to  be  proportional  to  [l-(Me/qMn)]  , where  Me  is  the 
entanglement  molar  mass,  Mn  is  the  number  average  molar  mass  of  the  polymer  before 
crazing,  and  q is  the  fraction  of  entangled  strands  that  do  not  undergo  chain  scission  in 
forming  the  craze.  Their  calculated  results  agreed  with  the  experimental  results  in 
polymethylmethacrylate  (PMMA)  from  Doll  (Doll,  1983). 

Chang  and  Hsu  (1991)  studied  the  effect  of  molar  mass  on  precrack  hysteresis 
energy  in  determining  the  ductile-to-brittle  transition  of  PC.  By  assuming  that  most 
inelastic  energy  is  used  for  creating  the  plastic  zone  ahead  of  the  crack  tip,  they 
established  the  relationship  between  hysteresis  loss  energy,  and  the  plastic  zone  volume. 
The  hysteresis  loss  energy  relates  in  almost  a linear  fashion  to  the  permanent 
displacement  of  the  samples.  This  plastic  zone  size  is  a function  of  deformation  rate, 
temperature,  specimen  thickness,  notch  radius,  deformation  displacement,  yield  stress, 
and  molar  mass.  The  results  show  that  higher  molar  mass  polymers  have  a higher 
precrack  hysteresis  energy  and  a larger  precrack  plastic  zone  volume  when  under  the 
condition  of  constant  loads. 

Effect  of  Viscoelasticity 

Viscoelasticity  is  a unique  property  of  polymers.  The  fracture  behavior  of 
polymers  is  complicated  due  to  this  property,  which  is  time-dependent.  The  response  of 
polymers  to  an  applied  stress  depends  upon  the  rate  of  loading  or  the  time  frame  of  the 
experiment.  If  the  rate  of  loading  is  high  or  the  time  frame  of  the  experiment  is  shorter 
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than  the  characteristic  relaxation  time  of  polymers,  they  will  display  more  elastic 
behavior.  On  the  contrary,  if  the  rate  of  loading  is  low  or  the  time  frame  of  the 
experiment  is  longer  than  the  characteristic  relaxation  time  of  polymers,  they  will  display 
more  viscous  behavior.  Hollmann  and  Hahn  (1989)  studied  the  effect  of  loading  rate  on 
the  fracture  toughness  of  epoxy.  Fracture  toughness  of  both  ductile  and  brittle  epoxies 
decreases  with  the  loading  rate.  Changing  the  test  temperature  also  alters  the 
characteristic  relaxation  time  of  polymers,  which  changes  their  response  to  an  applied 
stress  as  well. 

Effect  of  Type  of  Polymers 

The  fracture  mechanisms  of  polymers  are  different,  depending  on  the  types  of 
polymers.  The  fracture  mechanism  of  thermoplastics  can  be  the  breaking  of  the  primary 
covalent  bonds  or  the  secondary  Van  der  Waal  bonds.  If  amorphous  polymers  do  not 
have  enough  entanglement  or  semi-crystalline  polymers  do  not  have  enough  tie 
molecules,  then  their  fracture  mechanisms  will  be  the  chain  pulling  out  via  secondary 
bond  breaking.  This  mechanism  usually  requires  low  energy.  On  the  other  hand, 
polymers  with  high  entanglement  density  or  many  tie  molecules  will  fracture  via  primary 
bond  breaking,  which  requires  higher  energy.  However,  chain  scission  and  chain  pull  out 
may  occur  simultaneously  (Mohammadi  et  al.,  1992)  (Sambasivam  et  ah,  1995). 

Fracture  mechanism  of  thermosets  usually  occurs  via  primary  bond  breaking  due  to  the 
crosslink  points,  which  restrict  the  network  from  deforming  under  applied  stress. 
However,  energy  dissipation  in  polymer  fractures  is  complicated.  Nonelastic  behavior, 
hysteresis,  and  plastic  deformation  are  common  behaviors  in  polymer.  Shear  yielding 
and  crazing  are  two  competitive  micromechanisms  of  plastic  deformation.  Normally, 
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bulk  shear  yielding  will  result  in  ductile  fracture,  whereas  localized  shear  yielding  and 
crazing  will  result  in  brittle  fracture. 

Shear  Yielding  and  Crazing 

Shear  yielding  occurs  at  45  degrees  to  the  principal  stress  axis,  and  there  is  no 
change  in  volume.  Deformation  zones  by  shear  yielding  consist  of  polymer  drawn  to  a 
uniform  extension  ratio  with  no  fibrillated  structure  (Kramer,  1983).  Crazing  is  a 
cavitation  process  initiated  at  points  of  high  stress  concentrations.  Crazing  normally 
develops  in  a plane  perpendicular  to  the  maximum  principal  stress  and  causes  the  volume 
to  increase  in  the  system.  Crazes  consist  of  a network  of  highly  drawn  fibrils  spanning 
the  craze  surface,  separated  by  voids  (Plummer  and  Donald,  1989).  Typically,  crazing  is 
a precursor  of  a brittle  fracture,  since  it  is  highly  localized.  Energy  dissipation  from 
crazing  involves  a small  volume  of  materials.  However,  multiple  crazing  may  result  in 
ductile  fracture. 

Entanglement  density  has  significant  influence  on  the  competition  between  these 
two  mechanisms.  Kramer  (1983)  used  the  extra  surface  energy  hypothesis  to  compare 
PC  and  polytertbutylstyrene  (PTBS),  which  represent  the  two  extremes  of  entanglement 
density.  Assuming  that  only  chain  scission  is  involved  in  entanglement  loss,  the  surface 
energy  to  create  new  void  surface  area  can  be  estimated  by: 

T = y + 1 / 4dvEU... (2.1) 

where  T is  surface  energy,  y is  van  der  Waals  surface  energy  of  intermolecular  separation, 
d is  the  straight  line  distance  between  the  entanglement,  vE  is  the  entanglement  density, 
and  U is  the  energy  needed  to  break  a single  primary  chain.  By  letting  y and  U in  PC  and 
PTBS  to  be  the  same,  the  surface  energy  for  PC  will  be  higher  due  to  the  contribution 
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from  1/4 dvE.  Hence,  to  give  the  same  craze  tip  growth  rate,  the  craze  tip  stress  of  PC  will 
be  approximately  74  MPa,  while  the  craze  tip  stress  of  PTBS  is  27  MPa.  The  craze  tip 
stress  of  PC  is  higher  than  its  shear  yield  stress  at  room  temperature.  PC,  therefore,  tends 
to  deform  more  in  shear  yielding  than  in  crazing  at  room  temperature.  At  high 
temperatures,  the  rate  of  disentanglement  is  higher;  the  entanglement  that  is  lost  by  chain 
scission  is  smaller,  resulting  in  lower  surface  energy  and  also  in  craze  tip  stress.  Crazing 
becomes  more  competitive  for  PC  at  higher  temperatures.  On  the  other  hand,  with 
polymers  of  low  entanglement  density  such  as  PTBS,  in  which  entanglement  density 
contributes  little  to  the  surface  energy,  the  dependence  of  craze  tip  stress  may  be  less 
pronounced. 

There  are  two  types  of  polymers,  type  I and  type  II,  according  to  their 
micromechanism  of  plastic  deformation  (Fellers  and  Rahbar-Semanani,  1974).  Type  I 
polymers  contain  main  chain  aryl  units  connected  by  flexible  heteroatom  linkages,  such 
as  the  O-linkage.  Normally,  type  I polymers  are  ductile  in  their  amorphous  glassy  state. 
Examples  of  type  I polymers  are  bisphenol  A (PC)  and  polyethylene  terephthalate  (PET). 
Type  II  polymers  such  as  polystyrene  (PS)  and  PMMA  contain  rigid  units  in  the  main 
chain  and  no  flexible  linkages.  They  are  quite  brittle  in  their  amorphous  glassy  state. 
Plummer  and  Donald  (1989)  applied  the  concept  of  entanglement  to  differentiate 
between  these  two  types  of  polymers.  Type  I polymers  such  as  PC  or  polyethersulfone 
(PES)  have  a higher  entanglement  density  than  do  type  II  polymers  such  as  PS.  Surface 
craze  stress,  which  depends  on  the  entanglement  density,  is  higher  in  type  I polymers. 
Therefore,  type  I polymers  will  not  craze  at  room  temperature.  Transition  from  shear 
yielding  to  crazing  can  occur  when  the  temperature  increases.  This  transition 
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temperature  will  decrease  as  molar  mass  decreases.  Plummer  and  Donald  (1989) 
suggested  that  the  mechanism  of  crazing  in  type  I polymers  involves  the  motion  of  the 
whole  polymer  chains.  Their  data  also  showed  that  at  a lower  strain  rate  PC  and  PES  can 
craze  at  lower  temperatures  than  at  a high  strain  rate.  They  suggested  that,  as  the  time 
required  for  disentanglement  in  polymers  is  much  longer  than  the  time  required  for  chain 
scission,  disentanglement  plays  the  major  role  in  the  crazing  mechanism  of  type  I 
polymers. 

Wellinghoff  and  Baer  (1978)  studied  the  micromechanism  of  the  deformation 
process  of  type  I polymers,  type  II  polymers,  and  the  blend  of  them  as  well.  Deformation 
temperature  and  the  composition  of  the  blend  can  cause  the  transition  between  shear  and 
craze. 

Pitman  and  Ward  (1979)  calculated  the  craze  stress  of  PC  using  the  Dugdale 
Model  and  found  that  craze  stress  increases  as  a function  of  molar  mass.  Ishikawa  et  al. 
(1996)  studied  the  plastic  deformation  of  glassy  polymer  using  PC,  PMMA,  and  the 
blend  of  PC  and  PMMA  copolymer.  Their  results  showed  that  shear  yield  stress  is 
independent  of  molar  mass,  whereas  the  stress  for  nucleation  of  craze  increases  as  a 
function  of  molar  mass  in  PC  and  PMMA.  In  the  miscible  blend  of  PC  and  PMMA 
copolymer,  shear  yield  stress  and  the  stress  for  nucleation  of  craze  varies  as  a function  of 
the  content  of  PMMA  copolymer  following  the  composition  rule.  They  proposed  that  the 
mechanism  of  shear  yielding  involves  the  rotation  of  small  segments  in  polymer  chains, 
but  crazing  involves  a slippage  between  oriented  molecules  around  a void  as  a result  of 
high  dilatational  stress. 
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Fracture  Toughness  Concept 

Fracture  toughness  is  a measure  of  a material’s  resistance  to  crack  propagation. 
Unlike  tensile  strength,  which  is  flaw-dependent,  the  fracture  toughness  of  a material 
yields  a better  representation  of  material  characteristics.  There  are  two  approaches  to 
determine  fracture  toughness  of  a material:  the  energy  balance  approach  and  the  stress- 
intensity  factor  approach. 

Energy  Balance  Approach 

Griffith  (1921)  and  Orowan  (1948)  developed  the  concept  that  a fracture  occurs 
when  sufficient  energy  is  released  from  the  stress  field  by  growth  of  a crack  for  the  new 
fracture  surface.  There  will  be  a decrease  of  potential  energy  of  the  system  during  a 
fracture,  which  is  (7ta2a2t)/E,  according  to  Inglis  (1913)  in  the  case  of  an  infinitely  large 
plate  containing  an  elliptical  crack.  At  the  same  time,  there  is  an  increase  in  surface 
energy,  which  is  the  product  of  the  total  crack  surface  area  and  the  specific  surface 
energy. 

U -U0  = -( KG1a2t ) / E + 4atys...(2.2) 

where  U is  potential  energy  of  the  body  with  a crack,  U0  is  potential  energy  of  the  body 
without  the  crack,  cris  applied  stress,  a is  one-half  of  the  crack  length,  t is  the  thickness, 
E is  Young’s  modulus,  and  ys  is  specific  surface  energy  (energy  per  unit  area) 

By  rewriting  equation  2.2, 

U = 4 atys  - ( 7rcj2a2t ) / E + f/„...(2.3) 

At  equilibrium  condition, 


dU  / da  = 4 tys  - ( 2ncr2at ) IE-  0...(2.4) 
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2 ys  = {no1  a)  / £...(2.5) 

But  dU2 /da2  = -(2nd2 t)/E,  which  is  negative;  therefore,  the  equilibrium  condition 
described  by  equation  2.4  is  unstable  and  the  crack  will  always  propagate. 

By  rewriting  equation  2.5  for  plane  stress  (biaxial  stress)  condition, 

d = (2Eysl  naf'2... (2.6) 

This  condition  applies  to  a crack  in  a relatively  thin  plate,  in  which  the  stress  in  the 
thickness  direction  cannot  develop.  For  plane  strain  (triaxial  stress)  condition, 

a = [2Eys/na(\-v2)]U2...(2.7) 

where  vis  Poisson’s  ratio.  This  condition  applies  to  a crack  in  a thick  plate,  which  the 
stress  in  the  thickness  direction  can  develop.  Therefore,  a triaxial  stress  condition  is 
associated  with  the  suppression  of  strains  in  one  direction. 

Critical  strain  energy  release  rate,  or  Gic,  is  the  term  for  the  energy  required  to 
extend  a crack  over  a unit  area  for  mode  I or  tensile  mode  loading.  The  unit  of  Gic  is 
energy  per  unit  area  e.g.  J/m2.  At  equilibrium,  Gic  = 2ys\  therefore,  we  can  substitute  2ys 
in  equation  2.6  and  2.7  by  Gic- 
Stress-Intensity  Factor  Approach 

Based  on  LEFM,  which  applies  to  the  materials  that  obey  Hooke’s  law  or  that 
exhibit  a linear  relationship  in  stress-strain  curve,  and  by  methods  of  the  mathematical 
theory  of  elasticity,  stress  distribution  at  the  crack  tip  in  an  infinitely  large,  thin  elastic 
plate  with  a central  crack  can  be  derived  as  follows  (Mencik,  1992). 


(Tij(r,(p)  = a(a  / 2r)ul  fj(p...(2&) 
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where  a0(r,(p)  is  stress  distribution  at  the  crack  tip,  r,(p  are  the  coordinates  of  the  point 
considered,  cr  is  applied  tensile  stress,  a is  one-half  of  the  crack  length,  and/y(^)  is  a 
dimensionless  function  of  the  angle  between  the  crack  plane  and  the  respective  point. 
Figure  2.1  illustrates  this  stress  distribution. 

In  general,  there  are  three  basic  modes  of  loading  involving  different  crack 
surface  displacements  as  is  shown  in  Figure  2.2.  Mode  I is  the  opening  or  tensile  mode, 
where  the  crack  surfaces  move  directly  apart.  Mode  II  is  the  sliding  or  in-plane  shear 
mode,  where  the  crack  surfaces  slide  over  one  another  in  a direction  perpendicular  to  the 
leading  edge  of  the  crack.  Mode  III  is  the  tearing  or  antiplane  shear  mode,  where  the 
crack  surfaces  move  relative  to  one  another  and  parallel  to  the  leading  edge  of  the  crack. 
Mode  I is  the  most  simple  and  widely  used  method. 

Irwin  (1957)  determined  that  the  stress  field  around  a sharp  crack  in  a linear 
elastic  material  could  be  defined  by  Kj,  a stress-intensity  factor  which  is  a parameter  that 
is  independent  of  material  but  depends  on  the  magnitude  of  the  stress  field  and  the 
distance  from  the  crack  tip.  For  any  given  material,  fracture  will  occur  when  K]  exceeds 
Kic,  or  critical  stress-intensity  factor,  which  is  a material  property.  Subscript  I means 
tensile  opening  mode  or  mode  I,  and  subscript  C means  fracture  condition. 

Irwin  ( 1 962)  developed  the  relationship  between  stress-intensity  local  to  the  crack 
and  the  applied  stress  and  geometry  of  the  structure  during  loading. 


or 


Gii  = Kfy((p)/  (2;zr)1/2...(2.9) 
K = a(m)vl...(2.\0) 
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Figure  2.1  Stress  distribution  at  the  crack  tip  in  an  infinitely  large,  thin  elastic  plate  with  a 
central  crack.  Source:  Mencik  (1992) 
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Figure  2.2  Three  basic  modes  of  loading  and  crack  surface  displacements.  Source: 
Mencik  (1992) 
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Randall  (1967)  made  a review  of  Irwin’s  analysis,  which  applied  to  a flat  ellipse 
with  a sharp  crack  front  embedded  in  an  infinite  body  subjected  to  uniaxial  tension 
normal  to  the  plane  of  a crack,  as  illustrated  in  Figure  2.3,  and  showed  the  derivation  as 
follows: 

K = {<j / (j)\7ra)V2 {{a2  /c2)cos2  #>  + sin2  #>)1/4...(2.1 1) 
where  a is  half  of  minor  axis  of  the  ellipse,  c is  half  of  major  axis  of  the  ellipse,  and  ^is 
an  elliptical  integral  of  the  second  kind,  which  varies  between  1 for  a slit  crack  (a/b= 0)  to 
1.57  for  a/b=  1. 

At  (p  - tiJ2, 

K = a(m)V2  / </>...( 2.12) 

For  surface  cracks,  there  is  a correction  factor  of  front-face  free  surface  which  is 
about  a 20%  increase  in  K or  approximately  1 . 1 times  K.  As  a result,  for  surface  semi- 
elliptical cracks  of  depth  a and  half-width  c,  which  are  small  relative  to  the  thickness,  Ki 
(mode  I loading)  is  given  as  Mecholsky  (1993)  stated: 

K,  = l.lcrfl(;zx)1/2  / ^...(2.13) 

where  x is  a or  b,  whichever  is  smaller  (based  on  the  stress  distribution  at  the  end  of  the 
minor  axis  of  the  ellipse),  and  cra  is  applied  stress. 

Therefore,  at  fracture,  K/c  can  be  given  as 

KIC  = l.lcry (;zx)1/2  / ^...(2.14) 
where  <j/\s  the  stress  at  fracture. 

In  general,  Kj  may  be  defined  by 


K,  = Yaafn...(  2.15) 
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Figure  2.3  A flat  elliptical  crack  in  an  infinite  body  when  subjected  to  a uniaxial  stress 
normal  to  the  plane.  Source:  Randall  (1967) 
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and  KIC  = Yacf'n...(  2.16) 

where  Y is  a geometrical  constant  which  depends  on  shape,  location  of  flaw,  and  type  of 
loading,  cra  is  applied  stress,  crc  is  applied  stress  at  the  onset  of  crack  growth,  and/is  flaw 
size.  There  exist  tables  of  Y values  which  have  been  calculated  for  different  shapes  and 
locations  of  flaws  and  different  types  of  loading  (Murakami,  1987). 

Fractography  in  Polymers 

Fractography,  or  the  study  of  the  morphology  of  fracture  surfaces,  can  be  used  in 
post-failure  analysis  of  materials.  It  is  possible  to  identify  the  origin  of  a fracture,  its 
direction  and  speed  of  propagation,  the  fracture  stress,  and  also  the  mode  of  fracture  from 
studying  the  appearance  of  the  fracture  surfaces.  Research  has  been  done  on 
fractography  of  thermoplastics  (Mills,  1976),  (Matsumoto  and  Takemori,  1983),  (Chang 
and  Hsu,  1993),  (Chang  and  Hsu,  1994);  thermosets  (Yamini  and  Young,  1979), 

(Morgan,  1982),  (Kinloch  et  al.,  1983),  (Subramanian  et  ah,  1996);  and  the  blend  (Di 
Liello  et  ah,  1994)  using  different  tools,  e.g.,  the  optical  microscope  and  the  scanning 
electron  microscope.  The  different  features  of  fracture  surfaces  were  studied  in  order  to 
understand  the  various  fracture  mechanisms. 

Fractography  at  the  Macroscopic  and  Microscopic  Levels 

General  features  of  the  fracture  surfaces  of  brittle  fractures  are  shown  in  Figure 
2.4.  They  are  divided  into  4 regions. 

The  initiation  region,  or  defect  region,  is  the  region  where  the  fracture  has  started. 
The  fracture  origin  normally  is  due  to  machining,  impact,  or  material  defects  such  as 
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Figure  2.4  Typical  fracture  surface  of  brittle  materials.  Source:  Roulin-Moloney  (1989) 
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pores  or  inclusions.  The  material  defects  can  be  inherent  defects  that  exist  naturally  in  a 
material,  or  they  can  occur  as  a consequence  of  processing. 

The  mirror  region  is  an  area  with  a smooth  glossy  appearance  that  surrounds  the 
initiation  region.  The  crack  velocity  is  relatively  slow  in  the  mirror  region. 

The  mist  region  is  an  area  that  appears  to  be  substantially  rougher  than  the  mirror 
region.  The  crack  velocity  is  higher  than  in  the  mirror  region.  The  superposition  of  the 
primary  fracture  with  a secondary  fracture  can  create  parabolic  or  hyperbolic  markings  in 
the  mist  region.  The  open  ends  of  these  markings  indicate  the  direction  of  the  primary 
crack  front.  The  level  of  these  markings  is  different  from  the  surrounding  area.  The 
ridge-like  rims  have  been  proposed  to  be  the  result  of  viscous  flow  (Roulin-Moloney 
1989). 

The  hackle  region  appears  to  be  the  roughest  area  due  to  highest  crack  velocity. 
Macroscopic  crack  branching  initiates  at  the  end  of  the  hackle  region. 

In  the  case  of  ductile  fracture  surfaces,  they  are  also  composed  of  the  originating 
flaw,  but  the  mirror,  mist,  hackle  regions  are  not  as  obvious.  Normally,  the  cracks 
propagate  from  the  originating  flaw  at  a slow  speed,  then  the  speed  increases  until  crack 
branching  occurs.  Abrupt  transitions  can  be  seen  on  the  fracture  surfaces  between  the 
regions  of  slow  crack  speed  and  high  crack  speed  (Roulin-Moloney,  1989).  Rib 
markings  and  striations  are  commonly  found  on  fracture  surfaces.  Both  of  these 
markings  are  composed  of  crazed  materials.  The  former  is  due  to  multiple  crazes;  the 
latter  is  due  to  a single  craze.  There  are  markings  on  fracture  surfaces  showing  evidence 
of  plastic  deformation. 
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Lee  et  al.  (1993)  studied  the  fractography  of  injection-molded  PC/acrylonitrile- 
butadiene-styrene  (ABS)  blend.  They  found  that  the  fracture  surfaces  of  the  samples  in 
which  tensile  force  was  applied  parallel  to  the  injection  molding  direction  displayed  a 
herringbone  feature  and  possessed  greater  toughness.  On  the  other  hand,  the  fracture 
surfaces  of  the  samples  in  which  tensile  force  was  applied  perpendicular  to  the  injection 
molding  direction  displayed  a reverse  herringbone  feature  and  had  lower  toughness. 

They  proposed  that  the  different  position  of  the  secondary  crack  initiation,  which 
intersects  with  the  main  crack  from  these  two  fracture  surfaces,  contributes  to  the 
different  features  observed. 

Kim  et  al.  (1993)  studied  the  kinetics  of  the  process-zone  developed  during  stress 
relaxation  of  single-edged  notch  PC  samples.  This  process-zone  is  formed  and  developed 
surrounding  the  crack.  The  energy  required  for  the  growth  of  process-zone  is  much 
higher  than  the  surface  energy  associated  with  crack  formation;  therefore,  it  can  provide 
significant  resistance  to  crack  growth.  The  process-zone  was  found  to  consist  of  two 
intersecting  families  of  shear  bands  with  an  orientation  of  ± 59°  to  the  x-y  plane  with 
variable  density.  The  process-zone  length  grows  fast  at  the  initial  loading,  then  it 
increases  slowly  to  an  apparent  equilibrial  size,  depending  on  the  displacement. 

Kim  et  al.  (1994)  studied  the  effect  of  weathering  scale  and  rate  of  loading  on 
fracture  toughness  of  PC.  Compact  tension  samples  were  used.  Two  types  of  fractures 
were  found  in  the  study:  brittle  fractures  and  ductile  fractures.  In  the  load-displacement 
curve  of  a brittle  fracture,  the  load  rises  proportional  to  the  displacement  up  to  a critical 
value,  then  it  drops  to  zero.  In  a ductile  fracture,  a pop-in  phenomenon  occurs.  The  load 
rises  proportional  to  the  displacement  up  to  a critical  value,  then  it  drops  sharply  to  about 
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10%  of  the  maximum  value.  After  that,  there  is  a monotonic  increase  of  the  load  with  the 
displacement.  Slow  crack  growth  is  observed  during  this  stage.  A fracture  surface 
corresponding  to  the  brittle  fracture  can  be  divided  into  four  zones;  as  shown  in  Figure 
2.5,  zone  I is  a ductile  fracture  that  appears  ahead  of  the  crack,  zone  II  is  a brittle  fracture 
originating  from  the  crack  ahead  of  the  notch,  zone  III  is  a rougher  region  resulting  from 
multiple  cracks  or  crazes,  and  zone  IV  is  a smooth  region  of  dynamic  crack  growth.  A 
fracture  surface  of  a ductile  fracture  is  shown  in  Figure  2.6.  Zone  II  consists  of  the  crack 
arrest  line.  There  is  also  a large  plastic  deformation  zone  adjacent  to  the  brittle  fracture 
zone.  Tearing  marks  separate  these  two  zones. 

DeVries  and  Homberger  (1989)  studied  the  fracture  surfaces  of  PS  and  PC  using 
the  three-point  bend  test  and  the  impact  test.  They  concluded  that  the  untreated  samples 
have  flaws  generated  near  the  edge  of  the  samples  while  the  quenched  samples  (by  ice  or 
liquid  nitrogen)  have  flaws  generated  farther  from  the  edge.  Quenched  samples  from  the 
impact  test  are  rougher  and  exhibit  the  textured-feathered  shape  while  the  untreated 
samples  are  smooth  and  glassy.  The  size  and  shape  of  craze  in  PS  samples  are  also 
different.  The  untreated  samples  have  larger  craze  while  the  quenched  samples  have 
smaller  and  more  parabolic  craze.  However,  the  fatigue  fracture  samples  of  PC  exhibit 
the  opposite  fracture  surfaces  from  the  others.  The  untreated  and  annealed  samples 
consistently  exhibit  rough  and  coarse  surfaces,  while  the  quenched  samples  have  a large, 
smooth  area  near  the  notch.  They  explained  that  multiple  crazes  formed  under  a fatigue 
load  while  a single  craze  formed  under  an  impact  load. 

Chang  and  Hsu  (1994)  studied  ductile-to-brittle  transition  in  PC  as  a function  of 
test  temperature  and  found  differences  in  fracture  surfaces.  For  the  brittle  fracture 
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Figure  2.5  Schematic  diagram  of  a brittle  fracture  surface  of  a compact  tension 
polycarbonate  sample.  Source:  Kim  et  al.  (1994) 
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Figure  2.6  Schematic  diagram  of  a ductile  fracture  surface  of  a compact  tension 
polycarbonate  sample.  Source:  Kim  et  al.  (1994) 
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surfaces,  a crack  tip  blunting  zone  and  numbers  of  striating  lines  were  present,  but  there 
was  no  lateral  contraction.  For  the  ductile  fracture  surfaces,  the  distorted  and  irregular 
surfaces  were  observed  with  lateral  contraction.  For  the  semi-ductile  fracture  surfaces, 
various  features  were  observed  due  to  ductile-to-brittle  transition,  which  occurred  in  the 
whole  fracture  process. 

Atsuta  and  Turner  (1982)  compared  the  fracture  surfaces  of  polyethylene  glycol 
dimethacrylate  (PEGDM)  and  PEGDM  doped  with  1%  PMMA.  The  results  showed  that 
the  fracture  surfaces  of  undoped  PEGDM  mainly  consisted  of  brittle  features,  i.e.,  mirror, 
mist,  and  hackle  regions  with  river  markings.  The  mist  regions  were  few  and  limited  in 
area  and  included  linear  features,  which  were  formed  by  brittle  fractures  modified  by 
localized  plastic  deformation.  PEGDM  doped  with  1 % of  PMMA  also  showed  similar 
fracture  surfaces,  but  it  displayed  more  extensive  linear  features  and  also  a wide  variety 
of  morphologies  within  these  linear  features.  They  suggested  that  the  breaking  of  the 
main  crack  front  to  microcracks  and  the  joining  of  these  microcracks  in  order  to 
propagate  further  create  these  linear  features.  They  also  suggested  that  these  features  are 
created  via  localized  plastic  deformation,  which  incidentally  occurred  more  in  PEGDM 
doped  with  1 % of  PMMA. 

Fractography  at  the  Molecular  Level 

Kulawansa  et  al.  (1992)  studied  fracture  surfaces  of  PMMA,  PS,  and  PC  using 
the  scanning  tunneling  microscope  (STM).  All  fracture  surfaces  were  coated  in  gold  to 
35-50  Angstroms  thick.  STM  observations  were  made  on  macroscopically  smooth 
portions  of  each  fracture  surface.  PMMA  that  failed  under  tension  at  room  temperature 
consists  of  fibril-like  structures,  which  are  craze  fibrils;  bumps,  which  are  formed  as  the 
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cracks  are  deflected  around  regions  of  oriented  polymer  at  the  head  of  each  fibril;  and 
ridges,  which  are  formed  by  rows  of  bumps  when  crack  momentum  obscures  the 
individual  bump  along  a row.  Deep  cavities  are  also  found  on  PMMA  fracture  surfaces. 
These  features  are  formed  when  the  crack  penetrates  the  craze  boundary  into  the 
underlying  craze.  Secondary  fractures  within  the  craze  occurred  and  produced  hyperbolic 
features  on  fracture  surfaces.  Steps  are  also  observed  on  PMMA  surfaces  in  irregular  and 
in  regular  patterns.  The  minimum  distance  between  steps  is  typically  about  ten  times  the 
distance  between  ridges,  which  are  the  smallest  observed  longitudinal  structures.  Similar 
scaling  is  also  observed  in  PC  and  PS.  At  77  K,  fracture  surfaces  of  PMMA  are  rather 
smooth  when  compared  to  those  observed  at  room  temperature.  The  regular  step 
structures  which  are  found  may  be  the  consequence  of  the  high  constraints  on  craze 
growth  due  to  low  temperature.  Fracture  surfaces  of  PMMA  with  double  torsion 
geometry  are  much  smoother  and  have  a higher  periodicity  than  fracture  surfaces  failed 
under  tension.  With  double  torsion  geometry,  crack  velocity  is  slower.  At  higher 
magnifications,  alternating  regions  of  very  smooth  and  rough  surfaces  are  observed. 

They  proposed  the  mechanism  of  discontinuous  crack  growth,  in  which  the  crack 
propagates  along  the  craze  boundary  and  then  arrests  while  crazes  still  thicken  and  grow. 
Cracks  then  jump  again  along  the  craze  boundary  and  propagate.  Fracture  surfaces  of  PS 
failed  under  tension  at  room  temperature  are  similar  to  PMMA  with  higher  vertical  relief 
and  well-defined  ridges  due  to  a higher  draw  ratio  of  PS.  PC  is  the  only  polymer  under 
investigation  that  exhibits  shear  deformation  macroscopically.  Fracture  surfaces  of  PC 
failed  under  tension  at  room  temperature  also  show  ridges,  bumps,  and  steps.  The  ridge 
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spacings  and  step  spacings  are  much  smaller  than  those  observed  in  PMMA  and  in  PS. 
Steps  can  be  a consequence  of  shear  deformation  during  loading  and  fracture. 

Atomic  Force  Microscopy  as  a Tool  for  Fractography 

In  1982,  the  first  scanning  tunneling  microscope  (STM)  was  built,  which  enabled 
scientists  to  observe  a conducting  and  semi-conducting  surface  on  an  atomic  scale.  The 
images  from  STM  are  converted  from  either  the  spatial  variation  of  the  tunneling  current 
or  the  spatial  variation  of  the  tip  height.  The  atomic  force  microscope  (AFM)  was 
invented  subsequently  in  order  to  detect  insulating  surfaces.  AFM  is  a powerful  tool  in 
studying  surface  properties  of  materials,  e.g.,  surface  roughness  and  fractal  dimension. 
AFM  has  two  modes:  the  contact  mode  and  the  TappingMode®.  In  the  contact  mode,  the 
principle  is  based  on  the  measured  repulsive  force  between  the  tip,  which  is  located  at  the 
end  of  a cantilever,  and  the  surface  of  the  sample  via  the  cantilever  deflection.  Generally, 
the  tip  is  made  of  silicon  nitride  in  different  shapes  and  sizes.  In  the  TappingMode®,  the 
force  is  kept  constant;  the  cantilever  deflection  corresponds  to  the  height.  The  force 
brought  to  the  sample  surfaces  in  the  TappingMode®  is  greatly  reduced  compared  to  that 
brought  in  the  contact  mode  due  to  the  short  time  that  the  tip  contacts  with  the  surface. 
The  normal  forces  are  in  the  order  of  fractions  of  a nanoNewton  and  shear  forces  are 
zero.  Therefore,  it  is  suitable  for  soft  surfaces  like  polymers. 

The  Concept  of  Fractals 

At  the  molecular  level,  fracture  surfaces  can  be  modeled  by  fractals  (Mandelbrot, 
1984).  "Fractal  is  a term  derived  from  the  Latin  adjective  fractus,  which  is  related  to  the 
verb  frangere,  or  to  break,  in  English"  (Mandelbrot,  1977,  p.  294).  It  is  defined  as  an 
irregular  pattern  that  cannot  be  described  by  Euclidean  geometry  (Mandelbrot,  1977; 
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1983).  Fractals  can  be  in  the  form  of  curves  or  surfaces.  The  examples  of  fractals 
observed  in  nature  are  coastlines,  snowflakes,  the  shapes  of  clouds,  and  the  shapes  of 
leaves.  Self-similarity,  scale  invariance,  and  discontinuity  are  characteristics  of  fractals. 
Self-similarity  means  that  multiple  features  on  the  surface  appear  the  same  (Mecholsky 
and  Freiman,  1991).  Scale  invariance  signifies  that  a feature  at  one  level  of 
magnification  is  related  to  another  feature  at  another  magnification  through  a scalar 
magnification  constant  (Mecholsky  and  Freiman,  1991).  Discontinuity  characteristic  is 
defined  in  the  same  way  as  in  non-Euclidean  geometry,  i.e.,  the  slope  cannot  be 
determined  from  any  points  in  the  curves  or  surfaces. 

According  to  Euclidean  geometry,  we  can  describe  the  length  as  a function  of  the 
area  or  volume  as  follows  (Birdi,  1993): 

L oc  AI/2  x Ei/3...(2.17) 

where  L is  the  length,  A is  the  area,  and  V is  the  volume. 

In  the  case  of  non-Euclidean  geometry  we  can  describe  the  similar  relationship  as 
follows  (Birdi,  1993): 

L oc  Ml/Z\..(2.18) 

where  L is  the  length.  Mis  the  physical  quantity,  and  D is  the  fractal  dimension.  The 
quantity  D can  be  fractional,  in  contrast  to  Euclidean  geometry  in  which  it  is  an  integer. 
Fractal  dimension  can  be  used  to  quantify  irregular  patterns.  The  value  of  fractal 
dimension  will  be  1 < D < 2 for  irregular  curves  or  a two-dimension  profile,  and  will  be  2 
< D < 3 for  irregular  surfaces  or  the  three-dimension  profiles  shown  in  Figure  2.7. 

There  are  several  methods  for  measuring  fractal  dimension,  e.g.,  the  box  counting 
method,  the  power  spectrum  method,  the  variation  method,  and  the  slit-island  method. 


D = 2.0 


Figure  2.7  Fractal  dimensions  of  Euclidean  and  non-Euclidean  geometry. 
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F ractal  dimension  can  be  determined  directly  by  measuring  dimension  quantity 
from  the  fracture  surfaces  or  by  using  image  analysis.  The  algorithms  used  in  this 
dissertation  are  the  variation  method  and  the  slit-island  method. 

Variation  Method 

Variation  method  is  based  on  the  concept  of  ^-oscillations,  which  is  the  localmax  - 
localmin  difference  on  the  surface  (Spanos  and  Irene,  1993).  By  varying  the  local  points 
(i,  j coordinates  from  the  AFM  image),  the  averaged  sum  of  these  ^-oscillations  for  all  N 
x N points  or  s-variation,  V/,  can  be  obtained.  This  expresses  the  volume  that  the  surface 
occupies  when  it  is  approximated  at  the  observation  scale  defined  by  k as  shown  in 
Figure  2.8.  The  procedure  of  calculating  ^-variation  is  repeated  using  different  tile  sizes 
or  observation  scales.  Fractal  dimension  is  the  slope  of  the  log-log  plot: 

log  [(1  / £ „ ) : 3 Vf  (£„ )]  = a [log(l  /£„)]  + M 2.19) 

N N 

where  Vf(en)»  (1  / N2)'Yj'Yjv f(i,j,sn)  , kn  is  tile  size  varied  from  1,  2,  3....kmax,  Vis  all 

> j 

data  points  of  the  surface,  e„  is  kJN,  a is  the  slope  of  the  plot,  and  b is  a constant. 
Slit-Island  Method 

Based  on  the  Richardson  equation,  which  describes  the  self-similarity 
characteristic,  fractal  dimension  (Dp)  can  be  defined  as 

L = Zo£-(O'-D)...(2.20) 

where  L is  the  measured  length,  area,  etc.,  E is  the  scale  of  the  measurement  or  the 
yardstick  length,  L0  is  a constant,  and  D is  the  topological  dimension.  For  a straight  line, 
Dp  = D = 1 , and  for  a smooth  surface,  Dp  = D = 2.  For  a curve  or  a rough  surface.  Dp  lies 
between  1 and  2 or  2 and  3,  respectively.  Thus,  as  E approaches  0,  L approaches  oo. 
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Figure  2.8  Schematic  diagram  of  variation  method.  Source:  Spanos  and  Irene  (1993) 
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Instead  of  measuring  the  area  of  the  fracture  surfaces  themselves  compared  to  the  scale  to 
determine  fractal  dimension,  Mandelbrot  et  al.  (1984)  created  the  following  experimental 
technique.  The  fracture  surface  was  plated  with  nickel  and  mounted  in  an  epoxy.  Then 
the  specimen  was  polished  parallel  to  the  plane  of  fracture.  The  ‘islands’  of  steel 
appeared  on  the  surface  surrounded  by  nickel.  By  assuming  that  the  coastlines  of  these 
islands  are  self-similar,  the  fractal  dimension  of  these  coastlines,  being  curves,  can 
represent  the  fractal  dimension  of  the  fracture  surfaces.  The  relationship  between  the 
area  and  perimeter  of  the  islands  were  found.  The  slope  of  the  log-log  plot  between  these 
two  quantities  yields  the  fractal  dimension.  The  incremental  of  fractal  dimension  (Dp  - 1) 
obtained  from  the  slit-island  method  should  be  the  same  as  the  incremental  of  fractal 
dimension  (Dp  - 2)  obtained  from  measuring  the  fracture  surfaces. 

Hilders  and  Pilo  (1997)  summarized  the  algorithm  of  the  slit-island  method  as 

follows: 

Dr  = [(/(log 24V  <i(log  2 P,)]  - 1...(2.21) 

1=1  1=1 

where  Df  is  the  fractal  dimensional  increment  and  A,  and  P,  are  the  area  and  the 
perimeter  of  the  ith  island  on  each  layer  containing  n islands. 

There  are  arguments  in  the  literature  as  to  how  reliable  the  slit  island  method  is  in 
determining  fractal  dimension.  Pande  et  al.  (1987)  stated  that  the  slit-island  technique  is 
fundamentally  flawed  due  to  the  discrepancy  in  the  value  of  fractal  dimension  from  the 
slit-island  method  and  profile  Fourier  analysis.  Meisel  (1991)  analyzed  the  perimeter- 
area  and  the  perimeter-yardstick  relation  carefully  and  concluded  that  the  slit-island 
method  is  still  reliable  in  the  appropriate  range  of  yardstick.  Any  comparison  should  be 
made  at  the  same  yardstick.  The  islands  should  have  self-similar  characteristics,  but  they 


32 


do  not  have  to  possess  the  same  form  factor  or  the  same  ratio  of  length  to  breadth.  They 
suggested  that  for  each  island  and  for  the  range  of  yardsticks,  the  root  mean  square 
(RMS)  values  and  the  deviation  of  log  perimeter  and  log  area  should  also  be  determined. 
Then,  least-squares  fitting  with  points  weighted  according  to  their  RMS  deviations  should 
be  performed  in  order  to  determine  more  accurate  fractal  dimensions. 

Imre  (1992)  analyzed  two  sets  of  islands,  one  with  a constant  ratio  of  the  sides, 
the  other  with  a varied  ratio  of  the  sides.  The  results  showed  that  a linear  perimeter-area 
relation  can  be  obtained  even  in  the  case  of  islands  with  various  shapes,  but  the  slope  of 
the  line  does  not  represent  the  fractal  dimension  of  the  surface.  Therefore,  the  slit-island 
method  should  be  reliable  only  with  geometrically  similar  islands. 

Instead  of  using  the  Richardson  equation,  Underwood  and  Banergi  (1986)  used 
another  fractal  equation,  which  is 

log  Rl  ( 3 7)  = log C - ( D - 1)  log  7. . . (2.22) 

where  C is  a constant,  RL  is  the  profile  roughness  parameter  (the  apparent  profile  length 
divided  by  its  projected  length),  D is  fractal  dimension,  and  7 is  a measuring  unit. 

Fracture  surfaces  of  AISI  4340  steel  were  vertically  sectioned  and  fracture  profile  plots 
were  obtained  at  different  measuring  units.  The  fractal  plot  between  log  Rl(u)  and  log  7 
deviates  from  linear  at  the  upper  and  lower  ends.  They  suggested  that  the  plot  exhibits 
reversed  sigmoidal  curve  (RSC)  with  different  slopes.  Conventional  D is  determined 
from  the  slope  from  the  central  part  of  the  fractal  plot.  DIP  is  the  slope  determined  from 
the  inflection  point  in  the  fractal  plot.  The  relationship  of  conventional  D and  Dip  with 
tempering  temperature  shows  the  same  trend  but  different  values.  Surface  roughness 
parameters  Rl,  Rs  (the  apparent  surface  area  divided  by  its  projected  area),  and  RP  (the 


33 


ratio  of  the  average  peak  height  to  the  average  peak  spacing)  also  were  determined.  All 
of  these  surface  roughness  parameters  and  fractal  dimension  show  similar  trends  with 
tempering  temperature.  Further,  they  linearized  the  RSC  fractal  plot  and  determined  the 
constant  slope  P for  the  profile  and  y for  the  surface.  By  modifying  the  Richardson 
equation,  they  derived  Dp  = 0+  1 and  Dr=  2y+  2.  They  also  found  the  relationship 
between  Dip  and  [3  as  Dip  = 1.1 286/1°  02754.  Dp  and  Dy  decrease  with  tempering 
temperature  from  200°C  to  500°C,  then  they  increase  with  tempering  temperature  until 
700°C.  This  relationship  is  slightly  different  from  the  relationship  between  D,  Dip, 
surface  roughness  parameters,  and  tempering  temperature  mentioned  above.  They  did 
not  discuss  this  difference.  Therefore,  the  relationship  between  fractal  dimension  and 
toughness  cannot  be  stated  clearly. 

Quantitative  Fractography 

Quantitative  fractography  is  one  of  the  means  by  which  information  from  fracture 
surfaces  can  be  obtained  quantitatively  and  related  to  fracture  behavior  both  at  a 
macroscopic  and  a microscopic  level.  Fracture  toughness,  fractal  dimension,  and  residual 
stress  can  be  determined  using  quantitative  fractography. 

Conventional  techniques  for  fracture  toughness  determination  consume  much 
time  in  the  preparation  of  unique  specimen  dimensions  and  shapes  and  typically  require 
special  measuring  devices  such  as  a strain  gauge  for  recording  precise  deflection.  Using 
quantitative  fractography,  fracture  toughness  (Kic)  can  be  determined  from  standard 
samples  failed  under  tension  or  flexure  by  combining  the  information  from  fracture 


surfaces  and  the  fracture  stress. 
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It  has  been  shown  by  Mecholsky  et  al.  (1977)  that  for  most  glass  structures,  the 
failure-initiating  crack  can  be  modeled  as  an  equivalent  semicircular  crack  of  radius  c 
where  c = ( ab)u 2.  By  this  assumption  and  by  letting  (j)=  1 .57,  the  equation  for  Kjc  can  be 
written  as 

KIC  = \. 24a  f{(ab)m}m... (2.23) 

which  will  be  applied  to  the  analysis  of  the  brittle  PC  under  investigation. 

In  some  cases,  plastic  deformation  can  occur  at  the  crack  tip.  Assuming  that  the 
radius  of  the  plastic  zone  is  small  compared  to  the  crack  size,  LEFM  can  be  applied  to  the 
system  with  some  modifications.  Randall  (1967)  reviewed  Irwin's  analysis  and  showed 
that 

(KIC)2  = {\.2x(af)2c}/[<?>2  -{0.21 2(a  f)2  !(a  ys)2)]...(2.24) 
where  ays  is  yield  stress,  which  will  be  applied  to  the  analysis  of  the  PC  with  ductile 
fracture. 

Mecholsky  (1991)  has  assessed  quantitative  ffactography  and  has  shown  the  three 
important  equations  for  numerical  analysis  on  the  fracture  surfaces  of  many  inorganic 
brittle  materials,  which  are  also  composed  of  the  same  four  main  regions. 

crc'/2  = constant.. .(2.25) 
a(rJ)U2  = Mj/Y=KBj...(  2.26) 
c / Vj  - constant.. .(2.27) 

where  rj  is  the  radius  from  the  flaw  origin  to  each  boundary  (ri  = mirror  region,  r2  = mist 
region),  cr  is  fracture  stress,  Mj  is  a constant  for  each  boundary,  c is  critical  crack  or  flaw 
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size,  Y is  a geometrical  constant  for  crack  shape,  and  KBj  is  stress  intensity  at  each 
boundary. 

Mecholsky  (1993)  found  that  in  a soda-lime-silica  glass,  Mj  and  KBj  by  analogy 
are  constant  regardless  of  whether  or  not  the  flaw  is  in  the  plane  perpendicular  to  the 
maximum  stress.  It  was  shown  that  the  KiC  value  is  independent  of  the  depth-to-width 
ratio  of  the  flaw  and  can  be  calculated  from  the  mirror  constant  or  stress  intensity  at  each 
boundary.  The  ratio  of  the  square  root  of  the  depth  and  width  of  the  flaw  to  the  mirror 
radius  is  constant  for  semi-elliptical  flaws.  Mecholsky  et  al.  (1974)  calculated  fracture 
energies  from  the  mirror  data  for  most  glasses  over  a wide  range  of  compositions.  The 
results  agreed  well  with  measured  values.  Conway  and  Mecholsky  (1989)  used  crack 
branching  data  measured  from  the  fracture  surfaces  to  determine  the  amount  of  residual 
stress  of  tempered  and  annealed  aluminosilicate  glass  and  soda-lime  glass.  The  results 
agreed  well  with  the  residual  stress  determined  by  birefringence. 

Not  much  research  has  been  focused  upon  applying  the  technique  of  quantitative 
fractography  to  polymers.  Reed  and  Walsh  (1994)  studied  tensile  properties  of  epoxy, 
cyanate-ester,  vinyl-ester,  and  polyester  resin  systems  at  different  temperatures  and 
applied  the  technique  of  quantitative  fractography.  Their  analysis  of  linear  regression 
plots  of  tensile  strength  and  the  square  root  of  flaw  size  showed  that  at  low  temperature, 
the  sensitivity  of  tensile  strength  to  flaw  size  was  twice  the  sensitivity  at  room 
temperature.  Plangsangmas  et  al.  (1999)  applied  this  technique  to  determine  fracture 
toughness  (K/c)  of  an  amine-cured  epoxy  and  an  anhydride-cured  epoxy  cured  at  different 
periods  of  time.  They  also  correlated  Kjc  to  the  molar  mass  between  crosslink  (Mc). 
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As  the  fractal  dimension  becomes  a parameter  used  to  quantify  fracture  surfaces, 
attempts  have  been  made  to  study  the  relationship  between  fractal  dimension  and  other 
material  properties  such  as  fracture  toughness  and  impact  strength.  The  goal  is  to  predict 
the  material  properties  or  the  fracture  mechanism  from  fracture  surfaces. 

The  following  fundamental  relationship  between  the  fractal  dimension  and  the 
fracture  toughness  of  a material  has  been  derived  by  Mecholsky  et  al.  (1989)  using 
alumina  and  glass  ceramics. 

KIC  = Ea0'nD*u 2 = Y(0)afc'/2...(  2.28) 

where  Kjc  is  the  critical  stress  intensity  factor  or  fracture  toughness,  E is  Young’s 
modulus,  a0  is  the  characteristic  length  of  fracture  plane,  D*  is  an  increment  of  fractal 
dimension,  Y(q)  is  a geometrical  constant  dependent  on  the  geometry  of  the  crack  and 
loading  condition,  qyis  the  stress  at  fracture,  and  c is  the  crack  size  at  fracture. 

Mecholsky  and  Freiman  (1991)  also  found  that  D*  = c/rj  where  rj  is  the  mirror  radius  on 
the  fracture  surface.  Mecholsky  and  Plaia  (1992)  found  that  a0,  which  is  calculated  from 
the  equation  2.28,  can  be  used  to  characterize  the  fracture  process  of  single-crystal  Si, 
polycrystal  ceramics,  glass,  and  glass  ceramics.  Thompson  et  al.  (1995)  used  Lithia- 
based  glass  ceramics  with  microcracks  and  non-microcracks  in  the  study  and  found  the 
same  relationship  between  Ai/Cand  D*.  D*  is  the  same  in  the  mirror,  mist,  and  hackle 
regions,  and  a0  is  the  same  for  both  categories,  i.e.,  they  have  the  same  fracture  process. 

Research  in  the  field  of  metals  produced  different  results  from  that  in  ceramic 
materials.  Mandelbrot  et  al.  (1984)  found  that  fractal  dimension  has  a negative 
relationship  with  room  temperature  impact  energy  when  using  maraging  steel  with 
different  thermal  treatments.  Fractal  dimension  value  from  the  slit-island  method  agreed 
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well  with  the  value  from  the  fracture  profile  analysis  based  on  Fourier  analysis.  Mu  and 
Lung  (1988)  found  that  fractal  dimension  has  a negative  relationship  with  Kjc  using  steel 
with  different  heat  treatments.  The  following  equation  was  developed 

In  KIC  = constant  +(\-DF)\nel  / 2.. .(2.29) 

where  KJC  is  the  fracture  toughness,  DF  is  the  fractal  dimension,  and  e,  is  the  yardstick 
length,  which  is  fixed.  Pande  et  al.  (1987)  measured  the  fractal  dimension  of  titanium 
alloy  with  various  amounts  of  zirconium  using  the  slit-island  method  and  the  vertical 
section  method.  They  correlated  fractal  dimension  with  dynamic  tear  energy,  which  is 
found  to  decrease  with  increasing  quantities  of  zirconium.  Fractal  dimension  determined 
from  both  methods  shows  the  same  trend,  i.e.,  the  negative  relationship  with  dynamic  tear 
energy.  However,  the  values  are  different.  Ray  and  Mandal  (1992)  used  the  slit-island 
method  to  determine  fractal  dimension  of  high-strength,  low-alloy  (HSLA)  steel.  They 
found  a positive  relationship  between  impact  energy  at  different  subzero  temperatures 
and  fractal  dimension.  Hilders  and  Pilo  (1997)  used  the  slit-island  method  to  determine 
fractal  dimension  of  a commercial  medium-carbon,  ferrite-pearlite  steel  with  different 
thermomechanical  treatments  and  correlated  it  to  impact  toughness  tested  at  room 
temperature.  The  results  show  a positive  relationship.  They  explained  the  discrepancy  in 
the  relationship  from  that  found  by  Mandelbrot  et  al.  (1984)  and  Ray  and  Mandal  (1992) 
in  their  studies  in  terms  of  different  fracture  micromechanisms.  The  first  stage  of 
fracture,  cleavage,  is  the  main  fracture  mode.  The  microvoid  coalescence  is  very  small. 
Df  increases  with  toughness.  The  second  stage,  microvoid  coalescence,  is  dominant. 

The  surface  roughness  was  minimal  even  when  the  energy  absorption  is  high.  Thus, 
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Df  decreases  with  toughness.  Huang  et  al.  (1990)  determined  the  fractal  dimension  of 
impact-fractured  surfaces  in  a CK45  steel  with  different  thermal  treatments.  They 
compared  the  results  obtained  by  the  slit-island  method,  the  vertical  section  method,  and 
the  secondary  electron  line  scanning  method.  In  brittle  materials,  in  which  plastic 
deformation  is  very  small,  the  ‘islands’  on  a horizontal  section  of  the  fracture  surface 
should  exactly  match  with  the  ‘lakes’  on  the  same  section  of  a matching  fracture  surface. 
Both  islands  and  lakes  can  be  used  as  measured  objects  to  determine  fractal  dimension. 
On  the  contrary,  ductile  materials  with  much  plastic  deformation  associated  with  the 
coalescence  of  voids  can  create  dimples  on  the  fracture  surfaces.  As  the  toughness 
increases,  the  size  and  depth  of  dimples  increase.  If  lakes  or  dimples  are  used  as 
measured  objects,  the  perimeter-area  ratio  will  decrease  with  toughness  or  the  fractal 
dimension  will  decrease  with  toughness.  If  islands  are  used  as  measured  objects,  fractal 
dimension  will  increase  with  toughness.  This  may  explain  the  discrepancy  in  the 
literature  regarding  the  relationship  between  fractal  dimension  and  toughness  using  the 
slit-island  method.  The  relationship  between  fractal  dimension  and  toughness  shows  the 
same  trend  among  the  slit-  island  method  using  islands  as  measured  objects,  the  vertical 
section  method,  and  the  secondary  electron  line  scanning  method;  i.e.,  fractal  dimension 
increases  with  toughness.  However,  the  fractal  dimensions  obtained  from  these  three 
methods  have  different  values  because  they  describe  the  roughness  of  fractured  surface 
from  different  angles.  Alexander  (1990)  analyzed  fracture  surfaces  of  three  aluminum 
alloys.  Each  one  had  two  different  heat  treatments.  Underaging  and  overaging  were 
utilized  in  order  to  obtain  different  fracture  toughness  and  fracture  surfaces  while  keeping 
yield  strengths  similar.  The  analysis  of  digitized  profiles  of  sections  of  fracture  surfaces 
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indicated  that  the  profiles  are  fractal  in  nature.  However,  different  sections  from  the 
same  sample  may  have  different  fractal  plots.  Fractal  dimension  was  determined  using 
the  average  step  size  method  and  the  divider  method.  The  results  from  both  methods 
were  similar,  with  the  average  step  size  method  giving  more  accurate  data.  No  apparent 
relationship  between  fracture  toughness  and  fractal  dimension  can  be  found  from  the 
study,  whereas  a definite  relationship  between  fractal  dimension  and  linear  roughness  can 
be  seen.  Fracture  surfaces  with  a high  fractal  dimension  have  high  linear  roughness.  Su 
et  al.  (1995)  studied  impact- fractured  surfaces  of  steel  normalizing  at  different 
temperatures  by  using  the  slit-island  method.  They  found  that  fractal  dimension 
increases  with  impact  toughness,  but  it  decreases  with  normalizing  temperature.  They 
applied  the  slit-island  concept  to  the  microstructure  of  steel,  which  contains  ferrite  and 
pearlite.  The  isolated  pearlites  were  assumed  to  be  ‘islands’  and  ferrites  were  assumed  to 
be  Takes.’  Fractal  dimension  of  pearlite  increased  with  pearlite  size  and  also  with 
normalizing  temperature.  A similar  trend  was  found  for  fractal  dimension  of  pearlite 
determined  from  the  Richardson  equation.  The  interesting  point  from  this  study  is  that 
pearlite  microstructure  is  fractal  and  the  fractal  dimension  of  pearlite  can  be  used  to 
characterize  grain  size.  The  fractal  dimension  of  pearlite  has  an  inverse  relationship  to 
fractal  dimension  of  fracture  surfaces.  At  high  normalizing  temperatures,  pearlite  size  is 
large,  and  their  fractal  dimension  is  high;  however,  fractal  dimension  of  surfaces  is  low, 
and  impact  toughness  is  low.  The  fundamental  natures  of  these  relationships  are  not  well 
understood. 

In  the  field  of  polymers,  limited  research  has  been  done  on  fractal  dimension. 
Zhenyi  et  al.  (1991)  studied  epoxy  and  the  composite  of  epoxy-alumina  particles  and 
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found  that  K/c  values  are  significantly  different  while  fractal  dimensions  are  not.  Chen 
and  Runt  (1989)  found  that  fractal  dimension  has  a positive  relationship  with  KIC  when 
using  PS.  Lyu  et  al.  (1994)  used  the  microcrack  shear  band  theory  to  calculate  fractal 
dimension  for  ductile  polymers  and  found  an  inverse  relationship  with  KIC.  These  results 
from  the  literature  are  not  comparable,  since  they  used  different  polymers  and  different 
methods  to  determine  fractal  dimension.  However,  there  is  a trend  that  the  relationship 
between  fracture  toughness  and  fractal  dimension  may  be  based  on  fracture  behavior,  i.e., 
brittle  fracture  or  ductile  fracture  rather  than  on  materials.  Zhang  et  al.  (1996)  and  Zhang 
et  al.  (1997)  measured  fractal  dimension  of  wear  particles  of  polyetheretherketone 
(PEEK)  using  the  Richardson  equation  and  the  slit-island  method,  respectively.  Fractal 
dimension  values  from  both  methods  remain  almost  constant  in  the  low  range  of  contact 
pressure,  then  they  increase  with  the  contact  pressure  in  the  high  range.  A similar  trend  is 
found  for  the  relationship  between  contact  pressure  and  time-related  wear  rate.  Hence, 
fractal  dimension  can  be  used  to  predict  the  wear  rate  and  also  to  study  the  wear 
mechanism  of  polymers. 


Physical  Aging 

Definition  and  Characteristics  of  Physical  Aging 

Amorphous  glassy  materials,  both  ceramic  and  polymeric,  are  used  extensively  in 
many  industries.  Their  applications  range  from  household  products  to  aerospace.  The 
level  of  manufacturing  technology  is  quite  advanced,  yet  a fundamental  understanding  of 
their  physical  characteristics  is  still  lacking.  Structural  relaxation  of  glassy  polymers  is 
one  of  the  important  characteristics  that  has  been  studied  for  decades  because  of  the 
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changes  in  some  of  their  properties  during  application  which  result  from  this  relaxation. 
Struik  (1978)  defined  isothermal  structural  relaxation  below  the  glass  transition  (Tg)  of 
materials  as  physical  aging.  These  two  terms  may  be  used  interchangeably.  Structural 
relaxation  is  a phenomenon  that  occurs  in  amorphous  glassy  materials  below  their  Tg. 
Upon  quenching  from  the  temperature  above  Tg,  the  materials  are  not  given  sufficient 
time  to  reach  the  equilibrium  state.  Consequently,  there  is  some  excess  thermodynamic 
quantity  such  as  volume/free  volume  and  enthalpy  in  the  structure  as  compared  to  the 
equilibrium  state.  If  the  glassy  materials  are  annealed  below  their  Tg,  there  will  be  some 
molecular  mobility  driving  the  structure  towards  equilibrium.  This  structural  relaxation 
is  a form  of  time-dependent  decay,  the  rate  of  which  decreases  gradually  with  time  at  the 
specified  temperature.  Therefore,  the  relaxation  time  in  this  sense  is  retardation  time. 
The  relaxation  of  materials  towards  equilibrium  manifests  itself  in  a change  in  the 
packing  structure  resulting  in  an  overall  densification.  Structural  relaxation  will  occur  in 
the  range  of  temperature  between  Tg  and  the  highest  secondary  transition  (Tp),  according 
to  Struik  (1978). 

Kovacs  et  al.  (1979)  determined  three  important  characteristics  of  structural 
relaxation,  which  are  nonlinearity,  asymmetry,  and  the  memory  effect.  Nonlinearity  of 
structural  relaxation  means  that  a different  initial  magnitude  of  deviation  from 
equilibrium  will  result  in  a different  time  to  relax  for  the  same  amount.  Structural 
relaxation  also  exhibits  asymmetry  with  respect  to  the  sign  of  AT  (the  difference  from 
equilibrium),  or  in  other  words,  it  is  path-dependent.  The  relaxation  time  will  be 
different  between  the  heating  and  cooling  paths.  The  cooling  path  results  in  faster 
relaxation  than  the  heating  path.  Memory  effect  means  that  the  past  thermal  history  of 
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glass  will  have  an  influence  on  present  behavior.  With  regard  to  these  characteristics,  an 
experiment  in  structural  relaxation  must  be  done  very  carefully  in  order  to  obtain  correct 
results. 

The  difference  between  physical  aging  and  chemical  aging  is  highly  critical. 
Physical  aging  is  thermoreversible  because  no  permanent  change  in  the  primary  bond 
structure  or  molar  mass  occurs.  If  the  temperature  is  raised  to  above  Tg,  there  will  be  a 
long-  range  order  of  molecular  mobility,  which  erases  all  of  the  physical  aging  effects. 
Chemical  aging  is  non-thermoreversible  because  of  a permanent  change  in  the  primary 
bond  structure.  Cooling  rates  for  most  polymer  processing  are  rapid;  therefore,  glassy 
polymers  in  which  the  service  condition  is  below  Tg  will  undergo  physical  aging.  Due  to 
this  phenomenon,  structural  changes  of  polymers  affect  properties  of  polymers  during 
various  applications. 

Mechanism  of  Physical  Aging 

The  mechanism  of  physical  aging  can  be  described  by  two  approaches:  the 
thermodynamic  approach  and  the  free  volume  approach.  In  the  thermodynamic 
approach,  as  the  glass  solidifies  from  supercooled  liquid,  the  rate  of  molecular  motion 
decreases  gradually  and  if  the  cooling  rate  is  high  enough,  there  is  a departure  of 
thermodynamic  properties  from  the  equilibrium  state.  Enthalpy  and  volume  are  the  two 
most  important  properties.  In  this  non-equilibrium  state,  free  energy  is  larger  than  it  is 
for  the  equilibrium  liquid  at  the  same  temperature,  pressure,  and  concentration;  therefore, 
there  is  a drive  of  these  excessive  thermodynamic  properties  towards  equilibrium.  Based 
on  this  approach,  structural  relaxation  can  be  monitored  via  enthalpy  relaxation  or 
volume  relaxation  which  occurs  during  annealing  below  Tg. 
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Enthalpy  relaxation  is  usually  studied  using  a calorimetric  method  such  as 
Differential  Scanning  Calorimetry  (DSC).  When  the  annealed  or  aged  polymer  is  heated, 
there  is  a recovery  of  enthalpy,  which  is  lost  during  the  structural  relaxation.  The 
endothermic  peak  at  Tg  is  distinct,  which  indicates  that  the  annealed  or  aged  polymer 
needs  more  energy  to  superpose  on  the  specific  heat  change  associated  with  glass 
transition  compared  to  the  quenched  one.  Change  in  enthalpy  is  defined  by  the  area  of 
the  endothermic  peak  of  the  aged  sample  compared  to  the  quenched  one.  Flick  and  Petrie 
(1980)  studied  the  aging  of  PC  at  120,129,136,  and  143°C  and  found  that  change  in 
enthalpy  increases  with  aging  time  until  reaching  an  asymptotic  value  at  each  annealing 
temperature.  Biddlestone  et  al.  (1991)  and  Kemmish  and  Hay  (1985)  studied  enthalpy 
relaxation  of  polyetherimide  (PEI)  and  amorphous  PEEK,  respectively.  The 
experimental  results  show  that  the  unaged  part  of  enthalpy  [( \-AHt)/AHmax ] decreases 
logarithmically  with  time  at  each  aging  temperature  whereas  the  relaxation  time  can  be 
determined  using  the  Arrhenius  plot. 

Besides  enthalpy  relaxation,  volume  relaxation  also  has  been  extensively  studied. 
Volume  relaxation  can  be  measured  directly  from  change  in  volume  by  using  a 
dilatometer  or  indirectly  from  density  measurement  via  Archimedes’  method  or  by  using 
a density  gradient  column.  Change  in  volume  also  depends  on  annealing  time  and 
temperature.  The  relationship  between  enthalpy  and  volume  relaxation  with  physical 
aging  of  different  amorphous  and  semi-crystalline  polymers  has  been  studied  extensively 
(Kovacs  et  al.,  1979;  Carfagna  et  al.,  1988;  Ogale  and  McCullough,  1987;  Aref-Azar  et 
al.,  1983;  Vigier  and  Tatibouet,  1993;  and  Chang  and  Brittain,  1982). 
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In  the  free  volume  approach,  according  to  free  volume  theory,  which  explains 
molecular  motion  in  the  bulk  state,  molecular  motion  takes  place  if  there  are  holes  or 
voids  in  the  system.  There  will  be  an  exchange  of  places  between  molecules  and  holes. 
In  terms  of  polymer  chains,  the  same  concept  can  be  applied,  but  more  than  one  hole  will 
be  required  according  to  the  cooperative  motion  of  the  segments  of  chains.  The 
redistribution  of  excess  free  volume  occurs  during  physical  aging.  The  rate  of  this 
distribution  is  determined  by  the  mobility  of  polymer  chains.  Meyer  et  al.  (1990)  studied 
free  volume  change  in  poly  vinylacetate  (PVAc)  subjected  to  physical  aging  using 
fluorescence  spectroscopy.  The  result  indicated  that  in  the  initial  stage  of  aging,  the 
fluorescence  emission  intensity  of  the  probe  Auramine  O dispersed  in  polymer  matrix  is 
proportional  to  the  characteristic  times  for  stress  relaxation  of  the  polymer.  But  in  the 
latter  stage  of  aging,  the  intensity  is  less  sensitive  to  changes  in  total  free  volume.  Heater 
and  Jones  (1991)  and  Hill  et  al.  (1990b)  studied  free  volume  change  in  PC  subjected  to 
physical  aging  using  Positron  Annihilation  Lifetime  Spectroscopy  (PALS).  They  found 
that  the  lifetime  of  the  orthoPositronium  (oPs),  which  indicates  open  volume  size,  and  the 
intensity  of  the  oPs,  which  indicates  the  concentration  or  the  amount  of  open  volume, 
both  decrease  after  annealing.  They  found  that  the  activation  energy  calculated  for  the 
aged  PC  is  greater  than  that  of  the  unaged  PC,  which  implies  that  there  should  be  an 
increase  in  hindrance  of  motion  of  the  groups  or  segments.  Sandreczki  et  al.  (1995) 
studied  physical  aging  of  PC  at  24°C  using  PALS  and  found  that  the  mean  free  volume 
hole  size  remains  constant,  but  the  amount  of  free  volume  holes  decreases  as  a function 
of  aging  time.  The  relationship  between  fractional  free  volume  holes  and  the  annealing 
time  can  be  described  by  the  Doolittle  model  which  is  used  for  volume  relaxation. 
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Washer  (1985)  studied  the  relationship  between  the  change  in  free  volume  and  the 
change  in  enthalpy  relaxation  of  PC  and  found  no  correlation  between  the  change  in 
volume  and  enthalpy  relaxation.  The  decrease  in  free  volume  is  insignificant  compared 
to  the  change  in  enthalpy  relaxation  or  mechanical  properties. 

The  argument  that  free  volume  change  due  to  physical  aging  is  the  effect  of 
structural  change,  rather  than  its  cause,  is  reliable.  The  real  mechanism  during  structural 
relaxation  may  be  based  on  conformation  change  of  groups  or  segments,  which  result  in 
changes  in  free  volume.  Changes  in  free  volume  alone  are  not  enough  to  explain  the 
physical  aging  mechanism. 

Model  for  Physical  Aging 

Researchers  have  studied  and  established  mathematical  models  for  physical  aging 
in  order  to  predict  its  behavior  when  it  is  subjected  to  aging  conditions.  Kovacs  et  al. 
(1979)  introduced  a transparent  multiparameter  theory  to  explain  the  mechanism,  which 
added  sets  of  ordering  parameters  to  temperature  and  pressure  to  determine  the  state  of 
glass.  Each  of  these  parameters  possesses  one  retardation  time.  These  retardation  times 
depend  on  temperature,  pressure,  and  the  instantaneous  state  of  the  system.  The 
contribution  of  each  ordering  parameter  and  its  retardation  time  to  the  departure  originate 
the  structural  retardation  spectrum  for  the  system.  Adam  (1965)  introduced  the 
molecular  kinetic  theory,  which  explains  the  temperature  dependence  of  cooperative 
relaxation  in  glass-forming  liquids.  The  size  of  cooperative  region  is  temperature- 
dependent.  This  size  can  be  described  by  configurational  entropy  of  the  system. 
Matsuoka  and  Quan  (1991)  extended  Adam  and  Gibbs’  theory  to  the  non-equilibrium 
glassy  state.  They  introduced  the  model  for  intermolecular  cooperativity  in 
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conformational  relaxation.  The  individual  segments  cannot  relax  independently  of  their 
molecular  neighbors.  The  distribution  of  relaxation  time  broadens  in  the  non-equilibrium 
state  below  Tg  due  to  different  activation  energies  for  different  domain  sizes.  It  is  certain 
that  there  is  not  a single  relaxation  time  for  polymers  but  a spectrum,  which  consists  of 
sets  of  relaxation  time  from  different  segments.  Some  segments  relax  rapidly,  and  some 
relax  slowly. 

The  time-temperature  superposition  principle  states  that  time-dependent 
properties  at  one  temperature  can  be  superimposed  on  data  at  another  temperature  by 
shifting  the  data  along  the  axis  of  log  time  scale.  Struik  (1978)  applied  this  principle 
using  time-aging  time  to  the  study  of  creep  response.  He  showed  that  the  whole  spectrum 
of  relaxation  time  changed  during  aging.  Aging  influenced  all  relaxation  times  by  the 
same  factor  due  to  the  superimposability  of  the  creep  curve.  He  also  defined  the  aging 
rate  as  /u,  which  is  the  double-logarithmic  shift  rate,  ju  is  equal  to  -d  log  a I d log  te  where 
a is  the  horizon  shift  factor  from  the  creep  test  and  te  is  aging  time,  /.i  is  approximately  1 
over  a wide  range  of  aging  time  for  various  polymers  on  the  condition  that  they  are  far 
from  equilibrium.  As  equilibrium  approaches,  /u  starts  to  decrease  and  will  be  zero  at  a 
temperature  close  to  or  at  Tg.  The  range  of  temperature  in  which  p is  approximately  1 
differs  for  each  polymer.  This  range  of  temperature  may  conform  to  the  range  of 
temperature  in  which  polymers  exhibit  ductility,  implying  that  both  aging  and  ductility 

need  mobility  of  the  system.  Aging  rate  may  be  defined  as  d ( G '/G  "o ) / d log  te,  where 

/ // 

GIG  o is  the  ratio  of  the  shear  modulus  at  aging  time  t to  that  at  specific  aging  time, 
which  resembles  t = 0,  and  te  is  aging  time.  Venditti  and  Gillham  (1992)  determined  the 
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aging  rate  of  PMMA  using  this  method  and  the  result  agrees  with  the  value  determined 
from  creep  response. 

Tool  (1946)  introduced  the  fictive  temperature  as  a parameter  to  describe 
structural  change  of  a system.  This  parameter  will  add  to  temperature  and  pressure  in 
order  to  determine  the  fundamental  nature  of  nonequilibrium  glass.  The  fictive 
temperature  (7/)  is  the  temperature  at  which  glass  is  in  an  equilibrium  state,  as  shown  in 
Figure  2.9.  As  the  liquid  above  Tg  is  instantaneously  quenched  to  77  below  Tg,  P}  will  be 
the  property  in  the  nonequilibrium  state  at  77  and  time  t0.  The  extrapolated  line  from  P, 
parallel  to  the  slope  of  the  curve  P vs.  T in  the  glassy  state  will  intersect  the  equilibrium 
line  at  Tf(t0).  If  the  glass  is  aged  at  77  for  the  time  tn,  the  fictive  temperature  will 
decrease  to  T/(t„).  The  derivative  of  the  fictive  temperature  with  respect  to  time  will 
correspond  to  sensitive  properties  that  relax  to  equilibrium.  As  stated  in  Tool’s  equation, 
dTJdt  = ( Tf  - T)/t(T,  T/)  where  ris  structural  relaxation  time  which  depends  on  both 
temperature  and  fictive  temperature.  Mijovic  (1994)  summarized  this  information, 
stating  that  there  are  three  groups  of  models  of  structural  relaxation  time  which  all  use  the 
fictive  temperature  introduced  by  Tool.  They  are  based  on  the  Arrhenius  equation,  free 
volume,  and  configurational  entropy,  respectively.  The  model  based  on  the  Arrhenius 
equation,  or  the  Narayaswamy  model,  is  widely  used  because  it  yields  the  results  which 
best  fit  to  experimental  data.  A model  based  on  the  Arrhenius  equation  is  shown  in 
equation  2.30. 

t - A exp[AA h */RT  + ( 1 - X)A h * / RT/]... (2.30) 
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Figure  2.9  Schematic  diagram  of  relaxation  of  any  property  (p)  of  glass-forming  liquid. 
Source:  Scherer  (1993) 
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where  r is  structural  relaxation  time,  X is  parameter  partitioning  the  activation  energy  into 
two  parts,  the  effect  of  temperature  and  of  structure,  respectively  (0<X<1),  A h*  is 
activation  energy,  7}  is  Active  temperature,  A is  a constant,  and  R is  a gas  constant. 

Effect  of  Physical  Aging  on  Mechanical  Properties 

Physical  aging  has  significant  effects  on  many  mechanical  properties  such  as 
creep  response,  yield  stress,  elongation,  impact  strength,  and  fracture  toughness. 
Generally,  creep  resistance,  yield  stress,  post-yield  stress  drop  (PYSD),  and  tensile 
strength  increase,  whereas  elongation,  impact  strength,  and  fracture  toughness  decrease 
with  aging  (Zurimendi  et  al.,  1982;  Adam  et  ah,  1975;  Golden  et  ah,  1967;  LeGrand, 

1969;  and  Hill  et  ah,  1990a).  Therefore,  it  is  important  to  take  the  physical  aging  effect 
into  account  in  design.  Struik  (1978)  has  studied  the  effect  of  physical  aging  on  creep 
response  of  various  polymers  and  used  the  horizontal  and  vertical  shift  factor  to  predict 
creep  response.  Some  work  has  been  done  in  order  to  study  the  relationship  between  the 
change  in  these  mechanical  properties  and  the  change  in  thermodynamic  properties. 

Flick  and  Petrie  (1980)  found  that  change  in  yield  stress  parallels  change  in  enthalpy  in 
the  study  of  PC. 

One  important  change  in  mechanical  properties  due  to  physical  aging  is  the 
ductile-to-brittle  transition,  which  has  been  studied  for  years  in  order  to  understand  the 
mechanism  and  to  predict  the  lifetime  of  products.  Jabarin  and  Loifgren  (1991)  found  that 
there  is  a direct  correlation  between  percent  of  samples  failing  in  brittle  behavior  and 
enthalpy  relaxation  in  physical  aging  of  PET  at  65°C.  An  abrupt  change  between  ductile 
and  brittle  failure  is  associated  with  a critical  value  of  the  enthalpy  relaxation.  Yang  et  al. 
(1996)  studied  the  ductile-to-brittle  transition  induced  by  physical  aging  of 
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po  lypheny leneoxide  film  and  found  that  both  aging  time  and  temperature  have  a strong 
effect  on  this  transition.  The  ductile-to-brittle  transition  may  be  defined  by  a sharp 
decrease  in  elongation,  impact  strength,  or  fracture  toughness.  Conformation  change  and 
the  competition  between  crazing  and  yielding  are  two  mechanisms  for  ductile-to-brittle 
transition  proposed  from  the  literature. 

During  structural  relaxation,  polymer  chains  rearrange  themselves,  resulting  in 
changes  in  packing  structure.  These  spatial  rearrangements  may  cause  changes  in 
conformation  of  groups  in  the  chain.  Bubeck  et  al.  (1984)  used  solid-state  variable 
temperature  nuclear  magnetic  resonance  (VTNMR)  in  studying  conformation  change  due 
to  physical  aging  of  PC  and  PEC.  The  results  show  that  in  a glassy  state,  the  resistance  to 
1 80°  phenyl  ring  flipping  decreases  as  temperature  increases,  but  the  resistance  to  1 80° 
phenyl  ring  flipping  increases  as  a function  of  physical  aging.  Fourier  transform  infrared 
spectroscopy  (FTIR)  measurement  shows  the  frequency  shift  of  carbonyl,  in-chain 
carbonate  vibration,  and  out-of-plane  ring  modes  but  not  in-plane  ring  modes.  The 
results  are  controversial  and  the  data  interpretation  is  not  well  established.  Hill  et  al. 
(1990b)  reported  the  increase  of  activation  energy  of  relaxation  as  a function  of  aging. 

The  value  of  the  activation  energy  conforms  to  the  secondary  relaxation  of  the  phenyl 
group  and  the  carbonate  group.  These  are  the  results  of  closer  packing  structure  of  the 
chain  due  to  physical  aging.  Stolarski  et  al.  (1994)  used  Raman  Difference  Spectroscopy 
to  study  the  conformation  change  of  PC  aged  at  1 30°C.  The  results  show  small 
frequency  shifts  of  several  absorbances  as  a function  of  aging  time.  The  sinusoidal  curve 
shows  the  transition  between  high-and-low  energy  conformations  for  specific  segments 
of  the  chain.  The  carbonyl  group  in  a low-energy  state  moves  to  a high-energy  state  to 
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accommodate  the  movement  of  the  high-energy  phenyl  ring.  Moore  and  O’Loane  (1981) 
studied  conformation  change  in  PET  aged  at  50°C  using  FTIR  and  found  that  the  ethylene 
glycol  group  changes  conformation  from  trans  to  gauche.  A decrease  in  intensity  of  out- 
of-plane  motion  bands  and  an  increase  in  intensity  of  in-plane  motion  bands  in  the 
aromatic  region  is  also  observed.  These  changes  become  distinct  near  the  onset  of 
embrittlement  of  PET.  There  should  be  a correlation  between  conformation  change  and 
embrittlement  due  to  aging.  To  study  conformation  change,  special  spectroscopy  is 
needed  and  the  change  is  based  on  the  particular  polymer. 

In  general,  ductile  polymers  exhibit  yielding  before  failure,  while  brittle  polymers 
usually  fail  without  yielding.  There  is  a certain  amount  of  plastic  deformation  that  occurs 
during  the  failure  mechanism  of  ductile  polymers.  However,  localized  plastic 
deformation  such  as  crazing  may  induce  brittle  failure  to  polymers.  Hill  et  al.  (1990a) 
studied  the  aging  of  PC  and  proposed  that  the  embrittlement  caused  by  aging  took  place 
from  the  competition  between  crazing  and  yielding.  The  as-extruded  yielding  stress  is 
lower  than  the  crazing  stress,  whereas  the  aged  yielding  stress  is  higher  than  the  crazing 
stress.  Aged  PC  fails  in  brittle  mode  or  crazing,  whereas  quenched  PC  fails  in  the  ductile 
mode  or  yielding.  Fracture  toughness  of  the  aged  PC  is  lower  than  in  the  as-extruded  ~ 

1 8.5%.  Fracture  surfaces  of  tensile  specimens  also  show  the  difference  between  ductile 
mode  of  quenched  PC  and  brittle  mode  of  aged  PC. 

Approaches  to  Control  Physical  Aging 

As  described  above,  physical  aging  is  an  important  phenomenon  of  which  a more 
fundamental  understanding  is  needed  in  order  to  control  its  effect  on  properties.  The 
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approaches  to  control  physical  aging  are  the  mechanical  approach  and  the  chemical 
approach. 

Mechanical  approach 

The  influence  of  mechanical  orientation  or  deformation  to  control  physical  aging 
has  been  studied.  Mukherjee  and  Jabarin  (1995)  studied  the  physical  aging  of  oriented 
PET  and  found  that  the  rates  and  extents  of  relaxation  decrease  with  increasing 
orientation.  PET  was  stretched  in  the  rubbery  state,  resulting  in  a decrease  of  amorphous 
content.  Therefore,  the  effect  of  aging  was  reduced.  Othmezouri-Decerf  (1994)  has 
performed  mechanical  de-aging  to  the  well-aged  specimen  at  room  temperature  and 
studied  the  influence  of  mechanical  de-aging  on  subsequent  physical  aging.  They  applied 
a small  deformation  to  the  well-aged  specimen  in  the  pre-yield  region,  then  let  the 
specimen  recover  or  age  at  room  temperature  for  different  periods  of  time.  Stress-strain 
behavior  and  the  mechanical  damping  spectra  of  these  specimens  were  measured.  The 
results  show  the  enhancement  of  two  peaks,  a and  a",  and  the  stress-strain  behavior 
resembles  the  quenched  one.  The  enhancement  of  a can  be  described  by  increasing  in 
structural  temperature  as  a result  of  mechanical  de-aging.  Both  the  effect  on  or  peak  and 
yield  strength  agree  with  the  model  established  by  Bauwens  (1986).  The  a” process  has 
not  yet  been  interpreted.  Its  amplitude  is  independent  of  the  structural  temperature  and 
can  only  be  affected  by  the  recovery  period  after  mechanical  de-aging. 

Chemical  approach 

The  modification  of  the  chemical  structure  of  the  main  chain  has  been  done  to 
control  physical  aging.  Bubeck  et  al.  (1984)  incorporated  ester  groups  in  PC  and  by 
varying  the  structure  and  amount  of  ester  groups,  the  resistance  to  embrittlement  due  to 
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physical  aging  was  controlled.  The  large  amount  of  terephthalate  ester  groups  improve 
the  resistance  to  embrittlement  compared  to  PC  or  PC  with  isophthalate  ester  groups. 
Increasing  terephthalate  ester  groups  increases  the  viscous  component  in  the  solid  state  in 
the  fa  regime,  resulting  in  a higher  fa  magnitude  and  a lower  PYSD.  LeGrand  (1996  ) 
studied  the  physical  aging  of  PC  end-capped  with  phenol,  paracumylphenol  and 
chroman-I.  He  found  that  both  molar  mass  and  the  bulky  end  group  affect  the  fracture 
stress,  especially  in  the  range  of  long-term  aging.  He  proposed  that  the  time  of 
embrittlement  depends  on  molecular  weight  and  annealing  temperature  and  also  the 
relaxation  time  at  infinite  temperature,  which  is  a function  of  activation  energy. 

The  Influence  of  End  Groups  on  Physical  Properties  and  Physical  Aging 

Limited  publications  are  available  on  the  effect  of  end  groups  on  mechanical 
properties  and  on  the  physical  aging  effect.  Askadskii  (1989)  studied  the  effect  of  end 
groups  on  the  properties  of  polymers  and  found  that  the  Van  der  Waals  volume,  the 
cohesion  energy,  and  the  molecular  refraction  of  different  end  groups  reach  the 
asymptotic  value  of  the  repeating  unit  value  when  n or  the  number  of  repeating  units,  is 
approximately  10-20.  However,  specific  types  of  end  groups,  such  as  bulky  end  groups 
or  small  groups  with  specific  intermolecular  interaction  such  as  an  OH  group,  may  affect 
properties  differently.  The  bulky  end  groups  affect  the  Van  der  Waals  volume  of  an 
averaged  unit,  the  cohesion  energy,  and  the  molecular  relaxation  since  they  occupy  more 
space  and  their  mobility  requires  higher  energy.  The  OH  group  has  a significant  effect  on 
the  cohesion  energy,  the  solubility  parameter,  and  the  surface  tension,  but  it  has  little 
effect  on  the  Van  der  Waals  volume  and  the  molecular  relaxation.  Mansour  et  al.  (1994) 
studied  the  dielectric  relaxation  of  polystyrene  with  side-chain  and  chain-end  labels  and 
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showed  that  both  chain-end  and  side-chain  are  strongly  coupled  with  the  segmental 
orientation.  They  also  showed  that  increases  in  the  lengths  of  end  groups  shifted  the 
relaxation  maxima  to  lower  frequencies.  Wilkinson  et  al.  (1991)  showed  that  polysulfone 
with  reactive  end  groups  enhanced  the  fracture  toughness  of  a bismaleimide/diallyl 
network  compared  to  that  found  with  unreactive  end  groups  due  to  good  interfacial 
adhesion  between  phase-separated  domains.  Marks  (1994)  showed  that  the  types  and 
concentrations  of  end  groups  in  bisphenol  A-tetrabromobisphenol  A copolycarbonates 
have  major  effects  on  the  copolymer  solution  washability  and  thermal  stability  due  to  the 
different  polarities  of  end  groups.  Fellers  et  al.  (1974)  studied  the  dynamic  mechanical 
properties  of  end  group-modified  PS.  The  results  show  that  the  intensity  of  the  /. 3 
transition  is  influenced  by  the  molar  mass  of  the  polymer,  and  also  by  the  size  and  the 
structure  of  end  groups.  The  J3  transition  intensity  decreases  with  the  molar  mass  of  the 
polymer.  At  the  same  range  of  molar  mass,  bulky  end  groups  enhance  the  /?  transition 
intensity  due  to  more  local  free  volume.  On  the  other  hand,  end  groups  with  a very  long 
group  attached  will  decrease  the  / 3 transition  intensity,  since  the  movement  of  the  chain  is 
restricted. 

The  Relationship  between  Secondary  Relaxation  and  Physical  Properties 

The  secondary  relaxations  are  phenomena  involving  molecular  mobility  on  a 
small  scale  compared  to  that  observed  in  the  glass  transition.  They  are  characteristic  for 
specific  molecular  groups.  However,  significant  changes  in  physical  properties  of 
amorphous  polymers  going  through  the  secondary  relaxations  are  observed.  Liu  and 
Donovan  (1996)  studied  the  effect  of  secondary  relaxations  on  the  fracture  toughness  of 
semicrystalline  nylon  6,  amorphous  nylon  6IcoT,  and  their  blends.  The  plots  between 
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fracture  toughness  and  temperature  of  the  three  materials  show  similar  changes  around 
the  secondary  relaxations.  Shuster  et  al.  (1994)  studied  the  antiplasticization  effect  of  PC 
and  PCL  (polycaprolactone)  blend.  At  20  % of  polycaprolactone,  the  blend  has  higher 
modulus  and  yield  strength  but  lower  elongation  compared  to  pure  PC.  The  secondary 
relaxation  of  the  blend  is  suppressed  but  the  activation  energy  of  the  secondary  relaxation 
is  slightly  higher  compared  to  pure  PC.  By  adding  PCL,  the  blend  has  more  efficient 
packing  than  pure  PC,  resulting  in  a decrease  of  free  volume.  The  density  of  the  blend 
also  increases  compared  to  pure  PC,  although  PCL  itself  has  a lower  density  than  pure 
PC.  Jho  and  Yee  (1991)  studied  the  secondary  relaxation  motion  in  PC  and  concluded 
that  the  secondary  relaxation  of  PC  is  due  to  the  cooperative  motion  of  several  repeating 
units  along  the  chain.  They  also  showed  that  this  molecular  motion  is  affected  by 
changes  in  the  environment.  Schartel  and  Wendorff  (1995)  studied  the  secondary 
relaxation  behavior  of  a set  of  polyarylates  and  the  model  systems  composed  of  a single 
repeat  unit  and  end  groups  using  dielectric  relaxation.  All  systems  under  investigation 
exhibited  similar  secondary  relaxation  with  almost  the  same  activation  energy.  The 
differences  are  the  relaxation  strength  and  the  relaxation  time  distribution.  They 
concluded  that  the  secondary  relaxation  of  polyarylates  involves  both  inter-  and 
intramolecular  motions.  Bartos  et  al.  (1996)  studied  the  free  volume  microstructure  of 
PC  and  tetramethylpolycarbonate  (TMPC).  The  results  show  that  even  TMPC  has  high 
free  volume,  but  it  is  more  brittle  than  PC  due  to  the  lack  of  internal  flexibility  of  the 
chain.  The  benzene  ring  flip  is  inhibited  by  the  methyl  groups.  The  ductility  correlates 
well  with  the  secondary  relaxation,  which  indicates  the  ability  of  the  materials  to  absorb 
or  dissipate  energy. 
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Assuming  that  the  conformation  of  polymer  is  random  coil,  there  is  the 
probability  that  chain  ends  will  be  adjacent  to  some  segments  of  the  chain  and  will 
thereby  affect  the  molecular  motion  of  those  regions,  resulting  in  changes  in  shape  and 
magnitude  of  secondary  relaxation  and  in  ductility  as  well. 

Summary 

Ductile-to-brittle  transition  in  polymers  can  be  induced  by  changes  in  molar  mass, 
the  rate  and  temperature  of  testing,  and  thermal  history.  Fractography  at  the  macroscopic 
level  exhibits  different  features  of  ductile  and  brittle  fractures.  As  shown  in  the  literature, 
little  research  has  been  done  on  the  study  of  fracture  surfaces  at  the  molecular  level  or 
using  quantitative  fractography  in  polymers.  The  fundamental  relationship  between 
fractal  dimension  and  material  properties  is  still  controversial.  At  present,  the 
relationship  between  fractal  dimension  and  fracture  toughness  has  been  developed  for 
glass  ceramic  materials  only.  The  relationship  between  fractal  dimension  and  fracture 
toughness  of  polymers,  which  can  exhibit  either  brittle  or  ductile  behavior,  may  depend 
on  the  type  of  fracture  process. 

Physical  aging,  a phenomenon  that  can  induce  ductile-to-brittle  transition  in 
amorphous  glassy  materials,  can  be  used  to  make  a polymer  fracture  in  a more  brittle 
mode  than  when  it  is  in  an  un-aged  condition.  Furthermore,  the  idea  of  endcapping 
polymer  chains  with  bulky  groups,  which  is  believed  to  control  physical  aging  without 
changing  the  desired  mechanical  properties,  can  be  studied  by  comparing  the  aging  rate 
and  the  degree  of  ductile-to-brittle  transition  that  occurs  in  polymers  with  different  end 
groups.  The  effect  of  end  groups  on  secondary  relaxations  can  also  be  investigated  to 
determine  their  influence  on  ductility. 
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Since  yield  stress  increases  due  to  physical  aging,  the  ability  to  develop  the  plastic 
zone  should  decrease,  resulting  in  a decrease  in  plastic  zone  size,  which  can  be  studied 
from  fracture  surfaces. 


CHAPTER  3 

MATERIALS  AND  METHODS 
Materials 

Materials  used  in  this  study  are  bisphenol  A polycarbonate  from  General  Electric 
(GE)  Corporate  Research  and  Development,  Schenectady,  New  York.  The  typical 
chemical  reaction  and  the  structure  of  the  repeating  unit  are  shown  in  Figure  3.1.  Two 
different  groups  are  endcapped  to  the  polymer.  One  is  a paracumylphenol  group;  the 
other  is  a chroman-I  group.  The  chemical  structure  of  both  end  groups  is  shown  in 
Figure  3.2.  Typically,  end  groups  used  in  commercial  polycarbonate  are  phenol.  The 
molar  mass  of  the  samples  determined  by  GE  is  shown  in  Table  3.1 


Table  3.1  Molar  masses  of  polycarbonate  with  different  end  groups  used  in  the  study 


Sample  code 

End  groups 

Molar  mass  (kg/mole) 

92185 

paracumylphenol 

26.3 

92183 

paracumylphenol 

31.2 

92178 

paracumylphenol 

43.9 

22152 

chroman-I 

24.0 

92128 

chroman-I 

31.7 

22141 

chroman-I 

42.7 

Polycarbonate  in  this  study  was  made  by  GE  by  using  interfacial  polymerization 
of  bisphenol  A and  phosgene.  The  typical  diagram  of  polycarbonate  interfacial 
polymerization  is  shown  in  Figure  3.3.  The  polymerization  took  place  at  the  interface 
between  the  organic  and  aqueous  phases.  Phosgene  was  used  to  connect  bisphenol  A to 
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CH3  0 


ch3 


Bisphenol  A Phosgene 


-HCl 


ch3 


ch3 


Polycarbonate 


Figure  3.1  Typical  chemical  reaction  of  bisphenol  A polycarbonate. 
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paracumylphenol 


\^X)H 


Chroman-I 


Figure  3.2  Chemical  structure  of  end  groups  used  in  the  study. 
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PC  pellets 


Figure  3.3  Typical  diagram  of  polycarbonate  interfacial  polymerization.  Source:  Schnell 
(1964) 
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itself  and  a base  was  used  to  capture  the  released  hydrochloric  acid  (HC1). 
Bischloroformate  of  the  alkyl  phenol  was  added  to  the  reaction  to  obtain  better 
endcapping.  The  two  phases  are  separated  and  a dilute  aqueous  HC1  solution  was  passed 
to  the  organic  phase  in  order  to  remove  the  base.  The  organic  phase  was  re-separated  and 
washed  several  times  with  distilled  de-ionized  water.  Recovery  of  the  polymer  from  the 
organic  solution  can  be  done  either  by  precipitation  using  methyl  alcohol  or  hexane,  or  by 
flashing  off  the  organic  solution  using  steam.  Polymers  were  injection  molded  to  dog- 
bone  shape  according  to  ASTM  D638-91  Type  I. 

Methods 

Gel  Permeation  Chromatography 

Gel  Permeation  Chromatography  (GPC)  was  used  to  determine  the  molar  mass  of 
the  samples  used  in  the  study.  The  instrument  consists  of  a Waters  600E  System 
controller,  a Waters  610  Fluid  Unit,  a Waters  717  Auto  Sampler,  a Waters  410 
Differential  Refractometer  (DRI),  and  a Waters  996  Photodiode  Array  (PDA)  Detector. 
The  column  system  consists  of  four  columns  of  crosslinked  polystyrene  (Phenogel  5) 
connected  in  series.  The  columns  are  300  mm  in  length  and  4.6  mm  in  diameter.  The 
pore  size  of  crosslinked  polystyrene  in  each  column  is  100  A,  500  A,  50,000  A,  and 
500,000  A,  respectively.  The  temperature  of  the  system  was  maintained  at  35°C.  The 
flow  rate,  the  injection  volume,  and  the  run  time  were  set  at  0.40  ml/min.,  25  pi,  and  45 
minutes,  respectively.  The  data  was  collected  on  a computer  and  analyzed  by  using  the 
Mellennium™  Chromatography  Manager  software  from  Millipore  Corporation,  Milford, 
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Samples  were  taken  randomly  from  dog-bone  samples  and  dissolved  in 
tetrahydrofuran  (THF)  in  4-ml  autosampler  vials.  The  approximate  weight  of  sample  was 
20  mg.  The  samples  were  filtered  using  0.45  pm  Whatman  TF  25  syringe  filters  and  10 
ml  disposable  syringes  at  least  twice  prior  to  the  analysis.  The  samples  were  loaded  into 
the  autosampler  tray  and  the  run  was  programmed. 

The  instrument  was  calibrated  by  eight  polystyrene  standards  with  narrow 
distribution  molar  masses  of  860,  410,  194,  87,  28,  10.2,  3.55,  and  1.35  kg/mole  and 
toluene  of  molecular  weight  92  g/mole.  The  polydispersity  index  of  these  eight 
polystyrene  standards  are  1.06,  1.06,  1.06,  1.06,  1.05,  1.07,  1.05,  and  1.07,  respectively. 
The  coefficient  of  determination  (R2)  of  the  third-order  regression  between  log  molar 
mass  and  retention  time  is  0.9998.  The  standard  error  of  the  regression  line  is  0.0539. 
Differential  Scanning  Calorimetry 

A Seiko  DSC  220C  interfaced  with  Seiko’s  SDM/5600H  computer  was  used  to 
measure  the  glass  transition  of  PC  samples.  The  instrument  was  calibrated  for 
temperature  and  enthalpy  using  Indium  (In),  Tin  (Sn)  and  Bismuth  (Bi).  The  heating  rate 
was  10°C/min  and  the  flow  rate  of  dry  nitrogen  was  100  ml/min. 

All  DSC  experiments  were  performed  by  using  the  open  type  aluminum  pan  with 
a crimped  flange  cover  at  the  heating  rate  of  10°C/min  and  the  flow  rate  of  dry  nitrogen 
of  100  ml/min.  The  approximate  weight  of  samples  was  10  mg.  Alumina  powder  was 
used  as  the  reference.  All  scans  were  baseline-corrected  using  an  empty  pan.  The  glass 
transition  temperature  was  measured  as  the  midpoint  of  the  intersection  of  the  thermal 
response  for  the  glassy  and  liquid  states.  The  average  glass  transition  temperature  was 
obtained  from  the  second  heat  and  the  third  heat  of  each  measurement. 


64 


Thermal  Treatment 

All  samples  were  dried  in  the  vacuum  oven  at  121°C  for  4 hours  and  cooled  to 
room  temperature  in  the  vacuum  overnight,  then  they  were  kept  in  the  desiccator  prior  to 
any  thermal  treatment.  The  dried  samples  were  equilibrated  in  a drying  oven  at  Tg+25°C 
for  40  minutes  in  order  to  remove  any  residual  physical  aging  effects,  and  then  they  were 
air  quenched  to  room  temperature.  These  are  designated  as  ‘quenched’  samples  in  this 
study.  Some  of  the  quenched  samples  were  aged  in  a drying  oven  at  Tg-20°C  for  seven 
days,  then  they  were  air  quenched  to  room  temperature.  These  are  designated  as  ‘aged’ 
samples  in  this  study.  All  samples  were  kept  in  the  desiccator  prior  to  the  test. 

Dynamic  Mechanical  Spectroscopy 

A Seiko  DMS  200  tension  mode  analyzer  interfaced  with  Seiko’s  SDM/5600H 
computer  was  used  to  determine  the  dynamic  mechanical  properties  of  materials 
described  in  this  study.  The  instrument  was  calibrated  using  the  maximum  value  of  the 
dynamic  loss  modulus  (E")  at  1 Hz  for  the  glass  transition  temperature  of  PMMA  which 
is  1 17°C.  The  elastic  response  of  the  instrument  was  calibrated  using  tensile  modulus  of 
PC  measured  at  an  elongation  rate  10%/min  in  uniaxial  tension.  The  approximate  width 
and  thickness  of  all  samples  were  4 mm  and  0.5  mm,  respectively.  The  gage  length  was 
10  mm  for  all  samples.  The  furnace  was  purged  with  dry  nitrogen  at  the  flow  rate  of 
200ml/min  during  testing. 

DMS  experiments  were  performed  at  the  heating  rate  of  0.75°C/min,  frequency  at 
0.1,  0.5,  1,  5,  and  10  Hz.  for  activation  energy  analysis.  The  heating  rate  of  2°C/min  and 
frequency  at  1 Hz.  were  used  to  determine  the  influence  of  end  groups  on  dynamic 
mechanical  properties.  The  samples  were  isothermally  aged  in  the  DMS  to  monitor 
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changes  in  E'  and  E"  during  aging  and  to  measure  aging  rates  as  a function  of  aging 
temperature. 

Density  Measurement 

The  densities  of  PC  samples  were  determined  by  the  displacement  method 
according  to  ASTM  D792-91.  The  liquid  used  for  displacement  was  distilled  water.  Five 
samples  of  each  category  of  samples  were  measured.  The  approximate  dimension  of 
samples  was  1 5 x 15x3  mm.  The  water  adsorption  by  the  PC  samples  was  determined 
to  be  negligible  by  water  adsorption  studies. 

Uniaxial  Tensile  Test  at  Room  Temperature 

An  Instron  Universal  Test  Instrument  was  used  to  load  and  fracture  the  samples  in 
uniaxial  tension.  The  room  temperature  stress-elongation  properties  of  the  samples  were 
measured  according  to  ASTM  D638M-91 . Sample  dimensions  corresponded  to  Type  I in 
ASTM  D638-91.  The  approximate  thickness  was  3 mm  and  the  grip-to-grip  distance  was 
3.9  in.  (90  mm)  for  all  samples.  The  crosshead  speed  was  0.5  in.  (12.7  mm)/min  for 
every  test.  Three  to  eight  samples  of  each  category  were  tested. 

Fractography 

Fracture  surface  analysis  by  optical  microscopy 

The  fracture  surface  of  each  tensile  sample  was  analyzed  using  an  Olympus 
optical  microscope  at  50x  and  lOOx  magnification.  The  depth  a and  half-width  b of 
fracture-initiating  flaws  and  the  radius  of  the  surrounding  characteristic  topography,  i.e., 
the  mirror  and  the  mist  regions  on  the  fracture  surfaces,  were  measured  using  a reticle. 

The  magnitude  of  the  length  was  calibrated  using  a standard  glass  scale  at  each 
magnification.  The  size  of  fracture-initiating  flaws  c then  were  calculated  by  c = ( ab)]l 2. 
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Fracture  surface  analysis  by  atomic  force  microscopy 

The  fracture  surface  of  each  tensile  sample  was  analyzed  at  room  temperature  in 
air  using  the  Nanoscope®  AFM  from  Digital  Instruments,  Inc.  The  fracture  surfaces 
were  cut  with  a low-speed  diamond  saw  to  ~ 1 mm  thick  and  glued  to  the  0.5  in.-diameter 
steel  sample  holders  (pucks).  The  TappingMode®  was  used  with  the  large  scanner, 
which  the  maximum  scan  size  is  125  x 125  pm  and  the  maximum  Z-range  is  5 pm,  for  all 
analysis.  Silicon  nitride  tips  with  a high  aspect  ratio  were  used.  The  scan  rate  was  ~ 1 
Hz.  The  surface  area  used  in  the  analysis  was  ~ 800  x 800  nm.  The  resolution  was  512 
data  points.  The  images  were  flattened  and  planefitted  by  using  the  Nanoscope®  III 
software  equipped  with  the  AFM.  The  flatten  command  calculates  a least  square  fitted 
polynomial  for  each  scan  line  of  data  in  the  image  and  subtracts  that  polynomial  from  the 
scan  line.  This  process  eliminates  image  bow,  horizontal  bands,  and  other  unwanted 
information  from  the  image.  The  planefit  command  calculates  a single  first  order 
polynomial  for  the  point,  which  is  the  average  value  of  the  columns  when  the  X-axis  is 
selected  or  the  average  value  of  the  rows  when  the  Y-axis  is  selected.  The  command  fits 
a polynomial  to  the  points  and  subtracts  this  polynomial  from  each  line  in  the  image. 

All  images  were  analyzed  to  determine  the  roughness  and  fractal  dimension  by 
using  the  Nanoscope®  III  software  equipped  with  the  AFM  and  the  software  developed 
at  the  University  of  Florida,  respectively.  The  roughness  parameter  used  in  this  study  is 
the  root  mean  square  (RMS)  value  and  the  mean  roughness  value.  The  RMS  value  is  the 
standard  deviation  of  the  Z values  within  the  given  area  and  is  calculated  as  shown: 
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where  Rq  is  the  RMS  value,  Zave  is  the  average  of  the  Z values  within  the  given  area,  Z,  is 
the  current  Z value,  and  N is  the  number  of  points  within  the  given  area.  The  mean 
roughness  value  is  the  mean  value  of  the  surface  relative  to  the  center  plane  and  is 
calculated  as  shown: 

LyL, 

Ra  =1  !LxLy  ]\\f(x,y)\dxdy...(32) 

0 0 

where  Ra  is  the  mean  roughness  value,  f[x,y)  is  the  surface  relative  to  the  center  plane, 
and  Lx  and  Ly  are  the  dimensions  of  the  surface. 

Fractal  dimension  was  determined  by  using  the  variable  method  and  the  slit-island 
method.  The  variable  method  algorithm  is  described  in  equation  (2.19).  The  software 
used  in  this  calculation  has  the  following  steps: 

1 . Pick  the  point  i,  j from  the  surface  which  has  an  N x N two-dimensional  array 
of  points  (N  = 5 12  for  the  AFM  image). 

2.  Choose  the  tile  size  k (in  this  software  k=  1 . . . .40). 

3.  Calculate  the  difference  at  that  point  (/,  j),  i.e.,  localmax  - localmin  (s  - 
oscillation),  as  shown  in  Figure  2.8. 

4.  Repeat  step  3 for  aWNxN  points  of  the  surface. 

N N 

5.  Calculate  v/as  the  sum  of  all  these  e - oscillations  by  ^ vf  (/',  j,  en ) 

i j 

6.  Repeat  step  1 to  5 with  different  tile  sizes. 

The  examples  of  data  and  plot  are  shown  in  Table  3.2  and  in  Figure  3.4.  The 
slope  of  the  plot  in  Figure  3.4  is  the  fractal  dimension  of  the  fracture  surface. 
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Table  3.2  Examples  of  data  used  in  calculating  fractal  dimension  by  the  variation  method 


Col.l 

Tile  size(k) 
(length  unit) 

Col.  2 
In  (l/en); 
£n  — kn/N 

Col.  3 

In  [(l/sn)3  * (1/N2)  *col.  4] 

Col.  4 

N N 

XZM'’7’0 

i j 

1 

6.238 

24.575 

9.20E+07 

2 

5.545 

23.157 

1.78E+08 

3 

5.140 

22.329 

2.63E+08 

4 

4.852 

21.742 

3.46E+08 

5 

4.629 

21.286 

4.29E+08 

6 

4.447 

20.912 

5.10E+08 

7 

4.292 

20.596 

5.90E+08 

8 

4.159 

20.321 

6.69E+08 

9 

4.041 

20.079 

7.47E+08 

10 

3.936 

19.861 

8.25E+08 

11 

3.840 

19.664 

9.01  E+08 

12 

3.753 

19.483 

9.76E+08 

13 

3.673 

19.316 

1.05E+09 

14 

3.600 

19.162 

1.12E+09 

15 

3.530 

19.018 

1.20E+09 

16 

3.466 

18.883 

1.27E+09 

17 

3.405 

18.756 

1.34E+09 

18 

3.345 

18.636 

1.41E+09 

19 

3.294 

18.523 

1.48E+09 

20 

3.243 

18.416 

1.55E+09 

21 

3.193 

18.313 

1.62E+09 

22 

3.148 

18.216 

1.69E+09 

23 

3.103 

18.123 

1.76E+09 

24 

3.060 

18.034 

1.83E+09 

25 

3.020 

17.948 

1.90E+09 

26 

2.980 

17.866 

1.97E+09 

27 

2.943 

17.787 

2.04E+09 

28 

2.906 

17.711 

2.11E+09 

29 

2.871 

17.637 

2.18E+09 

30 

2.837 

17.567 

2.25E+09 

31 

2.804 

17.499 

2.31E+09 

32 

2.773 

17.433 

2.38E+09 

33 

2.742 

17.369 

2.45E+09 

34 

2.712 

17.306 

2.52E+09 

35 

2.683 

17.245 

2.58E+09 

36 

2.655 

17.186 

2.65E+09 

37 

2.627 

17.127 

2.71E+09 

38 

2.601 

17.071 

2.78E+09 

39 

2.575 

17.015 

2.84E+09 

40 

2.550 

16.961 

2.90E+09 
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In  (1/sn) 


Figure  3.4  Plot  to  determine  fractal  dimension  from  the  variation  method. 
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The  slit-island  method  is  based  on  equation  (2.20).  The  software  sliced  the 
fracture  surface  into  the  contour  line.  Then  the  software  used  different  yardsticks  from  3 
to  40  to  measure  along  the  contour  line  and  determined  the  number  of  counts  at  each 
yardstick.  The  log-log  plot  between  the  yardstick  length  and  the  number  of  counts  yields 
a negative  slope  due  to  an  inverse  relationship  between  the  yardstick  length  and  the 
number  of  counts.  The  slope  is  the  fractal  dimension  of  the  contour  line,  which 
represents  the  fractal  dimension  of  the  fracture  surface  as  well.  The  examples  of  data  and 
plot  are  shown  in  Table  3.3  and  in  Figure  3.5. 
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Table  3.3  Examples  of  data  used  in  calculating  fractal  dimension  by  the  slit-island 
method 


Col.  1 
Yardstick 
(length  unit) 

Col.  2 

Number  of  counts 

Col.  3 

Log  (col.  1) 

Col.  4 
Log  (col.  2) 

3 

437.00 

0.478 

2.641 

4 

335.25 

0.602 

2.525 

5 

259.20 

0.699 

2.414 

6 

220.50 

0.778 

2.343 

7 

187.29 

0.845 

2.273 

8 

162.13 

0.903 

2.210 

9 

144.22 

0.954 

2.159 

10 

129.00 

1.000 

2.111 

11 

117.00 

1.041 

2.068 

12 

107.25 

1.079 

2.030 

13 

98.54 

1.114 

1.994 

14 

91.44 

1.146 

1.961 

15 

85.34 

1.176 

1.931 

16 

78.70 

1.204 

1.896 

17 

74.47 

1.230 

1.872 

18 

70.22 

1.255 

1.847 

19 

66.37 

1.279 

1.822 

20 

62.45 

1.301 

1.796 

21 

59.34 

1.322 

1.773 

22 

56.94 

1.342 

1.755 

23 

54.84 

1.362 

1.739 

24 

52.60 

1.380 

1.721 

25 

50.04 

1.398 

1.699 

26 

47.97 

1.415 

1.681 

27 

45.73 

1.431 

1.660 

28 

44.70 

1.447 

1.650 

29 

43.04 

1.462 

1.634 

30 

41.56 

1.477 

1.619 

31 

39.67 

1.491 

1.598 

32 

38.91 

1.505 

1.590 

33 

37.21 

1.519 

1.571 

34 

36.27 

1.532 

1.560 

35 

34.83 

1.544 

1.542 

36 

33.96 

1.556 

1.531 

37 

32.98 

1.568 

1.518 

38 

32.66 

1.580 

1.514 

39 

31.29 

1.591 

1.495 

40 

30.49 

1.602 

1.484 
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Figure  3.5  Plot  to  determine  fractal  dimension  from  the  slit-island  method. 


CHAPTER  4 

RESULTS  AND  DISCUSSION 


Molar  Mass  Determination 


In  order  to  verify  the  molar  mass  value  from  GE,  GPC  was  used  to  analyze  all  the 
samples  used  in  the  study.  The  results  are  shown  in  Table  4.1. 


Table  4.1  Molar  mass  of  polycarbonate  with  different  end  groups 


Sample 

code 

End  groups 

Mn 

(kg/mole) 

Mw 

(kg/mole) 

Mz 

(kg/mole) 

Polydispersity 

index 

92185 

paracumylphenol 

12.1 

26.4 

42.3 

2.18 

92183 

paracumylphenol 

15.2 

31.2 

49.4 

2.05 

92178 

paracumylphenol 

21.0 

43.6 

69.3 

2.08 

22152 

chroman-I 

12.0 

25.9 

41.8 

2.17 

92128 

chroman-I 

16.1 

33.2 

52.5 

2.06 

22141 

chroman-I 

20.3 

44.3 

70.3 

2.18 

The  molar  mass  values  from  GE  in  Table  3.1  are  comparable  to  Mw  values  in 
Table  4.1.  The  discrepancy  between  the  GE  reported  values  and  the  measured  values  in 
this  study  is  based  on  the  different  solvent  used  in  the  experiments.  GE  used  chloroform 
as  a solvent  while  our  lab  used  THF.  However,  the  trends  are  similar.  The  molar  mass 
values  from  Table  4.1  will  be  used  in  this  study.  By  using  the  number  average  molar 
mass  from  Table  4.1,  the  amount  of  end  group  in  each  type  of  samples  can  be  calculated 
from  the  molar  mass  of  paracumylphenol  and  chroman-I,  which  are  212  and  268  g/mole, 
respectively.  The  amount  of  end  groups  for  each  type  of  sample  is  shown  in  Table  4.2. 
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Table  4.2  Amount  of  end  groups  in  each  type  of  sample 


Sample  code 

End  groups 

Mn  (kg/mole) 

Weight  percent  of  end  group 

92185 

paracumylphenol 

12.1 

3.51 

92183 

paracumylphenol 

15.2 

2.78 

92178 

paracumylphenol 

21.0 

2.02 

22152 

chroman-I 

12.0 

4.50 

92128 

chroman-I 

16.1 

3.33 

22141 

chroman-I 

20.3 

2.64 

It  may  be  difficult  to  see  the  influence  of  end  groups  on  properties  since  the  mass 
fraction  of  end  groups  is  very  small  compared  to  the  polymer.  Normally,  the  purpose  of 
endcapping  polymers  is  to  control  the  molar  mass  or  to  create  special  chemical  reactions 
at  the  end  of  the  chain.  However,  some  particular  properties  may  be  sensitive  to  the  size 
of  end  groups,  e.  g.,  the  relaxation  property.  Therefore,  the  idea  of  endcapping  polymer 
chains  with  bulky  groups  in  order  to  control  the  effect  of  physical  aging  but  retain  the 
important  properties  of  polymers  may  be  possible.  As  shown  in  Table  4.2,  the  weight 
percent  of  chroman-I  is  higher  than  that  of  paracumylphenol  at  the  same  range  of  molar 
mass,  so  the  influence  of  end  groups  may  be  more  prominent. 

Differential  Scanning  Calorimetry  Results 

The  glass  transition  temperature  is  highly  critical  in  physical  aging  experiments. 
To  conduct  physical  aging  experiments,  the  samples  must  be  aged  at  the  temperature 
below  their  glass  transition.  The  results  of  the  glass  transition  are  shown  in  Table  4.3. 
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Table  4.3  The  glass  transition  of  polycarbonate  with  different  molar  masses  and  end 
groups 


Sample 

Code 

End  groups 

Mw 

(kg/mole) 

T„  measured 
(°C) 

Tg  reported  from  GE 
(°C) 

92185 

paracumylphenol 

26.4 

138.7 

140.0 

92183 

paracumylphenol 

31.2 

140.6 

141.9 

92178 

paracumylphenol 

43.6 

148.3 

149.6 

22152 

chroman-I 

25.9 

145.3 

145.6 

92128 

chroman-I 

33.2 

147.3 

150.0 

22141 

chroman-I 

44.3 

151.5 

151.1 

The  discrepancies  of  the  values  between  Tg  measured  and  Tg  reported  from  GE  are 
due  to  the  different  heating  rates  in  running  the  DSC.  GE  used  the  heating  rate  of 
20°C/min  while  we  used  the  heating  rate  of  10°C/min.  However,  both  values  of  T are 
similar.  The  glass  transition  temperature  is  a thermal  property  that  is  based  on  the  rate  of 
heating  and  cooling.  When  the  rate  is  high,  polymer  chains  do  not  have  much  time  to 
respond  resulting  in  a higher  Tg.  The  experiment  in  this  study  will  be  based  on  the  Tg 
measured  in  the  laboratory.  Tg  also  depends  on  the  molar  mass,  as  shown  in  both  end 
groups.  The  higher  the  molar  mass,  the  higher  the  Tg,  as  is  well  known  from  the  Fox- 
Flory  equation 


Tg=Tga>-K/(aR-aG)Mn...(  4.1) 

where  Tg  is  the  glass  transition  temperature  at  a molar  mass  Mn  , Tgx  is  the  glass 
transition  temperature  at  infinite  molar  mass,  K is  the  material  constant,  aR  is  the  thermal 
expansion  coefficient  in  the  rubbery  state,  and  «g  is  the  thermal  expansion  coefficient  in 


the  glassy  state. 
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Considering  the  polycarbonates  at  the  same  range  of  molar  mass,  those  with 
chroman-I  end  groups  have  higher  Tg  than  those  with  paracumylphenol  end  groups.  The 
bulky  end  groups  create  a local  free  volume  around  the  chain  ends  which  is  very  small 
compared  to  the  free  volume  created  from  the  bulky  side  groups  along  the  chains. 
Therefore,  Tg  does  not  decrease  with  the  size  of  end  groups  as  it  decreases  with  the  size  of 
side  groups.  It  will  be  more  difficult  for  polymer  chains  with  bulky  end  groups  to  start 
the  motion  at  the  Tg,  resulting  in  a higher  Tg.  The  other  reason  that  may  explain  why 
bulky  end  groups  result  in  higher  Tg  is  the  volume  that  they  occupy  in  the  system.  Bulky 
end  groups  themselves  occupy  more  volume  than  small  end  groups.  Hence,  polymers 
with  bulky  end  groups  have  less  free  volume  left  per  unit  volume,  resulting  in  higher  Tg. 

Since  end  groups  affect  the  Tg  of  polymers,  it  may  help  in  processing  some 
polymers  with  high  Tg  by  modifying  the  end  groups  that  can  lower  the  Tg  without  altering 
other  main  properties.  However,  there  is  some  interaction  between  the  molar  mass  and 
the  size  of  the  end  group  in  terms  of  the  influence  on  Tg.  Based  on  the  experimental 
design  analysis  using  The  Minitab™  Software,  the  effect  of  end  group  and  molar  mass  on 
the  glass  transition  can  be  studied  at  the  same  time.  By  letting  the  end  group  be  variable 
A,  the  molar  mass  be  variable  B,  and  the  glass  transition  be  y-response  and  by  using  a 
two-level,  full  factorial  design,  the  design  matrix  for  the  study  can  be  generated  as  shown 
in  Table  4.4.  The  matrix  is  based  on  a single  replicate  of  y-response. 

The  high  levels,  and  the  low  levels  of  the  variables  are  represented  by  the 
numbers  1 and  -1,  respectively.  The  high  level  of  variable  A means  the  large  end  group, 
which  is  chroman-I.  The  low  level  of  variable  A means  the  small  end  group,  which  is 
paracumylphenol.  The  high  and  low  levels  of  molar  mass  for  polycarbonate  with 


77 


Table  4.4  Design  matrix  for  glass  transition  study 


Run 

A 

B 

Y 

Effect 

1 

-1 

-1 

138.7 

4.9  A 

2 

1 

-1 

145.3 

7.9  B 

3 

-1 

1 

148.3 

-1.7  AB 

4 

1 

1 

151.5 

paracumylphenol  end  groups  mean  an  Mw  of  43.6  and  of  26.4  kg/mole,  respectively.  The 
high  and  low  levels  of  molar  mass  for  polycarbonate  with  chroman-I  end  group  mean  an 
Mw  of  44.3  and  of  25.9  kg/mole,  respectively.  The  effect  of  end  group  and  molar  mass  on 
the  glass  transition  have  the  same  direction,  i.e.,  increasing  the  size  of  the  end  group  or 
increasing  the  molar  mass  will  increase  the  glass  transition.  The  effect  of  increasing 
molar  mass  is  more  significant  than  the  effect  of  increasing  end  group  size  as  shown  by 
comparing  the  values  of  7.9  and  4.9.  Yet  there  is  some  interaction  between  end  group 
and  molar  mass,  and  the  effect  of  the  interaction  has  a negative  value.  This  means  that 
increasing  the  size  of  the  end  group  and  molar  mass  at  the  same  time  may  result  in 
lowering  the  glass  transition  of  the  polymer.  This  preliminary  study  may  be  a guide  for 
modifying  the  glass  transition  of  the  polymer  without  altering  the  main  properties.  The 
main  effect  plot  and  the  interaction  plot  are  shown  in  Figure  4.1. 

Another  observation  from  the  glass  transition  study  is  that  all  polycarbonate  in 
this  study  is  amorphous.  There  are  no  exothermic  peaks  of  crystallization  or  endothermic 
peaks  of  melting  detected  in  the  DSC  between  25°C  to  350°C.  The  DSC  plots  showing 
the  glass  transition  of  polycarbonate  with  paracumyl  and  chroman-I  end  groups  are 
shown  in  Figures  4.2  and  4.3,  respectively.  The  plots  are  offset  in  the  y-axis  in  order  to 


improve  legibility. 
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Figure  4.1  Main  effect  plot  and  interaction  plot  of  the  glass  transition  study. 


DSC  (mW) 
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Figure  4.2  DSC  curve  of  polycarbonate  with  paracumylphenol  end  group,  10°C/min 
heating  rate,  100  ml/min  N2  flow  rate. 


DSC  (mW) 
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Figure  4.3  DSC  curve  of  polycarbonate  with  chroman-I  end  group,  10°C/min  heating  rate, 
100  ml/min  N2  flow  rate. 
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Density  Results 

Density  measurement  can  be  used  to  determine  the  packing  efficiency  of  polymer 
chains,  i.e.,  the  higher  the  packing  efficiency,  the  higher  the  density.  Physical  aging  can 
increase  the  small  amount  of  density  less  than  0.5  percent  (Struik,  1978).  The  results  of 
density  measurement  are  shown  in  Tables  4.5  and  4.6. 


Table  4.5  Density  of  polycarbonate  with  different  end  groups  and  molar  masses  in  the 
quenched  state 


Sample  code 

End  groups 

Mw  (kg/mole) 

Density  at  23°C  (g/cm3) 

92185 

paracumylphenol 

26.4 

1.1910±0.001 1 

92178 

paracumylphenol 

43.6 

1.1934±0.0005 

22152 

chroman-I 

25.9 

1.1856±0.0005 

22141 

chroman-I 

44.3 

1.1876±0.0008 

Table  4.6  Density  of  polycarbonate  with  different  end  groups  and  molar  masses  in  the 
aged  state,  aging  at  Tg-20°C  for  7 days 


Sample  code 

End  groups 

Mw  (kg/mole) 

Density  at  23°C  (g/cm3) 

92185 

paracumylphenol 

26.4 

1.1934±0.0012 

92178 

paracumylphenol 

43.6 

1.1938±0.0007 

22152 

chroman-I 

25.9 

1.1897±0.0012 

22141 

chroman-I 

44.3 

1.1929±0.0005 

As  shown  in  Table  4.5  and  in  Table  4.6,  the  density  increases  as  expected  as  a 
function  of  molar  mass  in  both  types  of  polycarbonate.  The  density  of  polycarbonate 
with  paracumylphenol  is  higher  than  that  of  polycarbonate  with  chroman-I  at  the  same 
range  of  molar  masses  in  the  quenched  state.  Upon  quenching  from  above  Tg, 
polycarbonate  with  larger  end  groups  cannot  pack  well,  resulting  in  more  free  volume  in 
the  system.  Hence,  the  density  is  lower.  The  results  also  show  that  after  aging,  all 
samples  are  denser  due  to  the  relaxation  of  the  free  volume  or  the  rearrangement  of  the 
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polymer  chains.  Density  measurement  can  be  related  to  free  volume  in  the  system 
indirectly.  In  order  to  determine  the  real  amount  or  the  size  and  distribution  of  free 
volume,  Positron  Annihilation  Lifetime  Spectroscopy  (PALS)  may  be  used. 

Bubeck  et  al.  (1984)  showed  that  the  percent  increase  in  density  does  not  conform 
to  the  embrittlement  caused  by  aging  in  PC  and  in  polyestercarbonate  (PEC).  PEC 
densified  more  than  PC  after  aging  at  the  same  time  and  temperature,  but  it  was  still  more 
ductile  than  PC.  An  increase  in  density  may  not  be  the  only  factor  involved  in  physical 
aging.  The  results  in  this  study  show  that  polycarbonate  with  chroman-I  densified  more 
than  polycarbonate  with  paracumylphenol.  The  percent  increase  in  density  after  aging  is 
shown  in  Table  4.7.  The  numbers  may  indicate  that  during  aging,  there  is  more 
rearrangement  of  chains  in  polycarbonate  with  chroman-I  than  in  polycarbonate  with 
paracumylphenol.  However,  no  conclusion  can  be  made  about  the  degree  of 
embrittlement. 


Table  4.7  Percent  increase  in  density  of  polycarbonate  with  different  end  groups  and 
molar  mass  after  aging 


Sample  code 

End  groups 

Mw  (kg/mole) 

Percent  increase  in 
density 

92185 

paracumylphenol 

26.4 

0.20 

92178 

paracumylphenol 

43.6 

0.03 

22152 

chroman-I 

25.9 

0.35 

22141 

chroman-I 

44.3 

0.45 

In  order  to  compare  the  effect  of  end  group,  molar  mass,  and  thermal  treatment  on 
density  at  the  same  time,  the  Minitab™  Software  is  also  used  to  generate  the  design 
matrix  and  to  calculate  the  effect  of  three  variables.  By  letting  the  end  group  be  variable 
A,  the  molar  mass  be  variable  B,  thermal  history  be  variable  C,  and  the  density  be  y- 
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response  and  by  using  a two-level,  full  factorial  design,  the  design  matrix  for  the  study 
can  be  generated  as  shown  in  Table  4.8.  The  matrix  is  based  on  a single  replicate  of  y- 
response. 

Table  4.8  Design  matrix  for  density  study 


Run 

A. 

B 

C 

Y 

Effect 

1 

1 

-1 

-1 

1.1856 

-0.0040  A 

2 

1 

1 

-1 

1.1934 

0.0021  B 

3 

1 

-1 

1 

1.1934 

0.0031  C 

4 

1 

1 

1 

1.1940 

0.0006  AB 

5 

1 

-1 

-1 

1.1910 

0.0016  AC 

6 

1 

1 

-1 

1.1876 

- 0.0002  BC 

7 

1 

-1 

1 

1.1897 

0.0008  ABC 

8 

1 

1 

1 

1.1929 

The  high  levels,  and  the  low  levels  of  the  variables  are  represented  by  the 
numbers  1 and  -1,  respectively.  The  levels  of  variable  A and  variable  B are  the  same  as 
in  the  glass  transition  study.  The  high  level  of  variable  C means  the  aged  state;  the  low 
level  of  variable  C means  the  quenched  state.  The  main  effect  plot  and  the  interaction 
plot  were  shown  in  Figure  4.4.  The  main  effect  plot  is  straightforward.  A small  end 
group,  high  molar  mass,  and  aging  increase  density.  The  effect  of  the  end  group  is 
greatest,  and  the  effect  of  molar  mass  is  smallest.  However,  the  interaction  effect 
between  end  group  and  thermal  history  is  also  important  since  the  effects  of  end  groups 
are  different  between  the  quenched  state  and  the  aged  state. 

Dynamic  Mechanical  Spectroscopy  Results 

DMS  experimentation  was  used  to  study  dynamic  mechanical  behavior.  The 
effect  of  molar  mass  and  the  end  group  on  dynamic  mechanical  behavior,  activation 
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Figure  4.4  Main  effect  plot  and  interaction  plot  of  density  study. 
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energy,  and  aging  rate  can  be  studied.  Figures  4.5,  4.6,  4.7,  and  4.8  show  the  effect  of 
molar  mass  on  the  DMS  plots  of  polycarbonate  with  both  end  groups.  There  is  no 
significant  difference  in  storage  modulus,  loss  modulus,  and  tan  8 in  the  glassy  region  due 
to  the  effect  of  molar  mass.  The  temperature  of  the  maximum  loss  modulus  and  the 
maximum  tan  5,  which  can  be  used  to  indicate  Tg,  are  in  the  same  trend  as  Tg  determined 
from  DSC.  The  effect  of  molar  mass  can  be  seen  in  that  Tg  is  higher  for  the  higher  molar 
mass. 

There  is  no  significant  effect  of  end  groups  on  storage  modulus  in  the  glassy 
region,  but  there  is  a difference  in  loss  modulus  and  in  tan  5 magnitude  (Figures  4.9  and 
4.10).  The  p-relaxation  around  -100°C,  which  is  typical  for  polycarbonate,  appears  in 
both  end  groups  but  with  different  features.  The  P-relaxation  of  polycarbonate  with 
chroman-I  is  broader  than  polycarbonate  with  paracumylphenol  of  equivalent  molar 
mass.  In  the  case  of  polycarbonate,  there  are  numerous  studies  of  p relaxation  at  -1 00  °C. 
Schaefer  et  al.  (1985)  concluded  that  this  relaxation  is  associated  with  cooperative  ring 
flip.  The  process  involves  at  least  two  rings  in  the  chains  and  also  4-5  rings  of  the  other 
chains.  Their  conclusion  may  point  to  the  effect  of  end  groups  as  well.  The  bulky  groups 
will  restrict  the  ring-flip  process  to  some  extent.  Fellers  and  Rahbar-Semanani  (1974) 
studied  the  dynamic  mechanical  spectroscopy  of  end  group-modified  polystyrene  and 
found  that  molar  mass  and  the  structure  of  the  end  group  affect  the  intensity  of  tan  delta 
in  the  P region  (60-80°C).  PS  with  a carboxyl  end  group  and  high  molar  mass  displays 
less  intensity  of  tan  delta.  Bulky  end  groups  enhance  the  p intensity  due  to  local  free 
volume,  whereas  long-chain  end  groups  decrease  the  P intensity  due  to  restriction  of  the 
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-e—  paracumylphenol  end  group,  Mw  43.6  kg/mole 


Figure  4.5  Storage  and  loss  modulus  of  polycarbonate  with  paracumylphenol  end  group 
as  a function  of  molar  mass,  2°C/min  heating  rate,  1 Hz.,  200  ml/min  N2  flow  rate. 


Log  Tans 
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Figure  4.6  Tan  6 of  polycarbonate  with  paracumylphenol  end  group  as  a function  of 
molar  mass,  2°C/min  heating  rate,  1 Hz.,  200  ml/min  N2  flow  rate. 


88 


Temperature  (°C) 


Figure  4.7  Storage  and  loss  modulus  of  polycarbonate  with  chroman-I  end  group  as  a 
function  of  molar  mass,  2°C/min  heating  rate,  1 Hz.,  200  ml/min  N2  flow  rate. 
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Figure  4.8  Tan  5 of  polycarbonate  with  chroman-I  end  group  as  a function  of  molar 
mass,  2°C/min  heating  rate,  1 Hz.,  200  ml/min  N2  flow  rate. 
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Figure  4.9  Storage  and  loss  modulus  of  polycarbonate  as  a function  of  end  group, 
0.75°C/min  heating  rate,  1 Hz.,  200  ml/min  N2  flow  rate. 
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Figure  4.10  Tan  8 of  polycarbonate  as  a function  of  end  group,  0.75°C/min  heating  rate,  1 
Hz.,  200  ml/min  N2  flow  rate. 
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movement  of  the  chain.  The  p region  they  studied  is  the  region  close  to  the  glass 
transition  and  its  shape  is  that  of  a shoulder,  not  a peak.  The  frequency  used  was  rather 
high,  i.e.,  110  Hz.  They  concluded  that  this  p transition  is  associated  with  the  oscillation 
of  end  groups.  The  origin  of  P-relaxation  can  be  the  motion  of  the  local  group  in  the 
main  chain  or  in  the  side  group.  The  results  of  this  study  suggest  that  there  is  some 
cooperative  movement  and  that  there  is  an  interchain  effect  on  P-relaxation.  Therefore, 
different  end  groups  can  create  different  features.  Chroman-I  creates  more  local  free 
volume  and  facilitates  the  movement  of  the  segments  nearby,  resulting  in  broader  P- 
relaxation. 

The  activation  energy  of  a-relaxation  at  Tg  and  P-relaxation  can  be  compared 
between  the  two  end  groups  by  using  DMS  with  multifrequencies.  The  activation  energy 
is  calculated  by  using  the  Arrhenius  equation  as  shown  in  equation  4.2 

/ = zfexp(_£/*n...(4.2) 

where /is  frequency,  A is  a constant,  E is  the  activation  energy  in  J/mole,  R is  a molar  gas 
constant  which  is  8.31434  J/mole/K,  and  T is  the  temperature  in  degrees  Kelvin.  DMS 
experiments  were  carried  out  at  0.75°C/min  heating  rate  and  frequencies  of  0.1,  0.5,  1,  5, 
and  10  Hz.  The  activation  energy  was  calculated  from  the  slope  of  the  straight  line  plot 
between  log  f and  1000/T.  Three  samples  were  run  for  each  end  group.  The  examples  of 
the  activation  energy  plot  for  a-relaxation  and  P-relaxation  are  shown  in  Figures  4.1 1 and 
4.12.  The  mean  and  standard  deviation  of  activation  energy  for  both  end  groups  is  shown 


in  Table  4.9. 
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Table  4.9  Activation  energy  of  a-relaxation  and  p-relaxation  of  polycarbonate  with 
different  end  groups 


Sample 

Code 

End  groups 

Mw 

(kg/mole) 

Activation  energy  at 
a-relaxation  (kJ/mol) 

Activation  energy  at 
P-relaxation  (kJ/mol) 

92183 

paracumylphenol 

31.2 

770.6±2.5 

42.6±4.9 

92128 

chroman-I 

33.2 

843.2±43.7 

32.2±2.0 

The  difference  in  activation  energy  between  these  two  end  groups  may  be 
explained  in  that  at  Tg,  it  is  more  difficult  for  polymer  chains  with  chroman-I  to  pass  each 
other  due  to  the  larger  size.  Therefore,  the  activation  energy  is  higher.  On  the  other 
hand,  at  p-relaxation  it  is  easier  for  polymer  chains  with  chroman-I  to  move  locally  due  to 
their  higher  local  free  volume.  This  difference  in  P-relaxation  may  result  in  different 
energy  dissipation  of  polymers. 

DMS  also  can  be  used  to  monitor  physical  aging  effects  in  polymers.  The  effects 
of  physical  aging  are  generally  seen  in  the  range  of  temperature  between  p-relaxation  and 
Tg  (Struik,  1978).  At  temperatures  below  P-relaxation,  polymer  chains  exhibit  very  little 
movement,  and  at  temperatures  above  Tg,  the  effects  of  physical  aging  are  erased. 

Figures  4.13  and  4.14  show  the  comparison  of  tan  8 plot  between  quenched  and  aged 
samples.  There  is  some  difference  in  tan  8 plot  for  both  end  groups  in  the  range  of 
temperature  between  -50°C  and  120°C.  In  order  to  verify  that  the  difference  does  not 
come  from  sample  variation,  the  same  samples  were  used  to  monitor  the  change.  The 
samples  were  equilibrated  and  air  quenched  as  described  in  Chapter  3,  then  they  were 
aged  in  the  drying  oven  at  Tg-40°C  for  7 days.  The  samples  were  put  in  DMS  and  run 
from  -150°C  to  70°C  at  2°C/min  and  1 Hz.  The  same  samples  were  removed  from  the 
DMS  and  equilibrated  and  air  quenched  then  run  in  DMS  again.  Figures  4.15  and  4.16 


Log  frequency  (Hz.) 
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Figure  4.1 1 Activation  energy  plot  at  a-relaxation.  Polycarbonate  with  chroman-I  end 
groups,  Mw=33.2  kg/mole 
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Figure  4.12  Activation  energy  plot  at  p-relaxation.  Polycarbonate  with  chroman-I  end 
groups,  Mw=33.2  kg/mole 
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Figure  4.13  Comparison  of  tan  6 of  polycarbonate  with  paracumylphenol  end  groups, 
Mw=3 1 .2  kg/mole  between  the  quenched  and  aged  states. 
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Figure  4.14  Comparison  of  tan  5 of  polycarbonate  with  chroman-I  end  groups,  Mw=33.2 
kg/mole  between  the  quenched  and  aged  states. 
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Figure  4.15  Tan  8 plot  of  polycarbonate  with  paracumylphenol  end  groups,  Mw=31.2 
kg/mole  as  a function  of  aging  using  the  same  sample. 
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Figure  4.16  Tan  8 plot  of  polycarbonate  with  chroman-I  end  groups,  Mw=33.2  kg/mole  as 
a function  of  aging  using  the  same  sample. 
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show  that  there  is  indeed  a difference  in  tan  6 in  the  temperature  range  of  -50°C  to  70°C. 
The  samples  were  run  only  up  to  70°C,  to  avoid  any  deformation  or  damage  to  the 
samples  that  could  occur  as  the  temperature  approached  T . 

Bubeck  et  al.  (1984)  compared  the  magnitude  of  tan  8 at  80°C  of  PC  and  PEC. 

The  difference  in  magnitude  was  ~ 0.019.  They  concluded  that  the  higher  the  magnitude, 
the  higher  the  resistance  to  aging.  The  tan  8 magnitudes  at  80°C  for  both  end  groups  are 
not  significantly  different.  They  are  in  the  range  of  approximately  0.007-0.014.  The 
origin  of  this  relaxation  in  PC  is  suggested  to  be  the  cooperative  movement  of  the  main 
chain  of  about  2-3  repeating  units.  The  samples  in  this  research  do  not  have  any 
difference  in  the  repeating  units  and  the  effect  of  the  end  group  on  this  relaxation  may  not 
be  enough  to  show  the  difference. 

Another  way  to  probe  the  effect  of  aging  is  to  determine  the  rate  of  aging  by  using 
DMS  experimentation.  The  samples  were  equilibrated  and  quenched  as  described  in 
Chapter  3.  The  samples  were  put  in  the  DMS  chamber  and  heated  to  the  aging 
temperature  at  5°C/min  and  maintained  at  that  temperature  for  at  least  400  minutes.  The 
change  in  storage  modulus  (E')  was  plotted  as  a function  of  log  aging  time.  The  slope  of 
the  linear  regression  line  is  the  rate  of  aging.  Aging  temperatures  were  Tg-20,  Tg-40,  and 
Tg-60°C.  Three  samples  were  used  for  each  aging  temperature.  As  shown  in  the  results 
(Figure  4.17),  the  average  rates  of  aging  for  both  end  groups  are  not  significantly 
different  from  each  other  at  Tg-Ta  = 40,  60°C  due  to  the  overlapping  of  their  standard 
deviations.  Only  at  Tg-Ta  = 20°C  did  the  rates  of  aging  of  these  two  end  groups  differ. 
Polycarbonate  with  paracumylphenol  end  groups  has  a higher  rate  of  aging  than 
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Figure  4.17  Rate  of  aging  of  polycarbonate  with  different  end  groups  as  a function  of 
aging  temperature. 
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chroman-I.  Ricco  and  Smith  (1990)  used  the  same  method  to  determine  the  rate  of  aging 
of  thin  film  commercial  polycarbonate  but  included  a constant  strain  of  2.6  percent  during 
the  experiment.  These  values  ranged  from  2.6  percent  at  30°C  to  3.5  percent  at  1 10°C. 
The  values  in  this  study  are  in  the  same  range  even  though  there  was  no  constant  strain 
applied  to  the  samples.  Struik  (1978)  determined  the  rate  of  aging  of  quenched  samples 
of  many  polymers  by  using  the  shift  rate  from  creep  compliance  experiment.  The  rate  of 
aging  of  polycarbonate  is  constant  from  -50°C  to  85°C  then  starts  dropping  when  the 
temperature  is  lower  than  -50°C  or  higher  than  85°C.  The  rates  of  aging  for  both  end 
groups  in  this  study  also  decrease  at  Tg-Ta  = 20°C,  which  is  higher  than  85°C. 

Tensile  Properties 

Tensile  tests  were  performed  for  quenched  samples  at  room  temperature. 
Polycarbonate  with  low  molar  mass  exhibits  brittle  fracture  with  low  elongation  and  no 
yielding,  while  polycarbonate  with  high  molar  mass  exhibits  ductile  fracture  with  high 
elongation  and  yielding.  Typical  stress-elongation  plots  are  shown  in  Figures  4.18  and 
4.19.  The  results  of  tensile  tests  for  quenched  samples  are  shown  in  Tables  4.10  and 
4.11. 

As  shown  in  Tables  4.10  and  4.1 1,  the  ductile-to-brittle  transition  induced  by 
differences  in  molar  mass  is  obvious.  More  entanglement  enhances  the  mechanical 
properties  of  polymers.  On  the  other  hand,  there  is  no  significant  difference  in  tensile 
properties  between  the  two  end  groups  at  the  same  range  of  molar  mass.  The  weight 
percentages  of  end  groups  are  very  small  compared  to  the  whole  polymer  chains.  They 
contribute  little  to  mechanical  properties. 
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Table  4.10  Results  of  tensile  test  for  polycarbonate  with  paracumylphenol  end  groups  in 
the  quenched  state 


Mw 

(kg/mole) 

Stress 

at 

yield 

(MPa) 

Elongation 
at  yield 
(%) 

Stress  at 
break 
(MPa) 

Elongation 
at  break 
(%) 

Young’s 

modulus 

(GPa) 

PYSD 

(%) 

26.4 

NA 

NA 

25.96 

4.1 

1.08 

NA 

(3.18) 

(0.5) 

(0.08) 

43.6 

62.68 

7.6 

56.31 

140.9 

1.54 

23.6 

(1.19) 

(0.4) 

(2.06) 

(10.7) 

(0.16) 

(3.1) 

Values  in  parentheses  represent  standard  deviation. 


Table  4.1 1 Results  of  tensile  test  for  polycarbonate  with  chroman-I  end  groups  in  the 
quenched  state 


Mw  Stress 

(kg/mole)  at  yield 
(MPa) 

Elongation 
at  yield 
(%) 

Stress  at 
break 
(MPa) 

Elongation 
at  break 
(%) 

Young’s 

modulus 

(GPa) 

PYSD 

(%) 

25.9 

NA 

NA 

22.74 

3.7 

1.06 

NA 

(2.92) 

(0.6) 

(0.06) 

44.3 

60.71 

10.2 

59.15 

164.2 

1.67 

23.2 

(0.59) 

(2.4) 

(0.30) 

(3.2) 

(0.03) 

(1.2) 

Values  in  parentheses  represent  standard  deviation. 


Another  set  of  quenched  samples  were  aged  in  the  drying  oven  at  Tg-20°C  for  7 
days  in  order  to  study  the  effect  of  physical  aging  on  tensile  properties.  Tensile  tests 
were  also  performed  at  room  temperature.  The  results  of  tensile  tests  of  the  aged  samples 
are  shown  in  Tables  4.12  and  4.13. 

As  shown  in  the  results,  physical  aging  can  induce  ductile-to-brittle  transition  in 
the  samples  with  high  molar  mass.  Yield  stress  increases,  breaking  stress  decreases, 
elongation  decreases,  and  post-yield  stress  drop  (PYSD)  increases.  Physical  aging  affects 
the  packing  efficiency  of  polymer  chains,  which  relates  to  the  rearrangement  of  polymer 
chains.  Decrease  in  free  volume  is  a consequence  of  this  phenomenon.  Therefore, 
yielding  can  be  affected  dramatically.  Polymers  require  higher  stress  to  yield  due  to  their 
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Figure  4.18  Room  temperature  stress-elongation  plot  of  polycarbonate  with 
paracumylphenol  end  groups  in  the  quenched  state.  Crosshead  speed  is  0.5  inch/min. 
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Figure  4.19  Room  temperature  stress-elongation  plot  of  polycarbonate  with  chroman-I 
end  groups  in  the  quenched  state.  Crosshead  speed  is  0.5  inch/min. 
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Table  4.12  Results  of  tensile  test  for  polycarbonate  with  paracumylphenol  end  groups  in 
the  aged  state 


Mw  Stress 

(kg/mole)  at  yield 
(MPa) 

Elongation 
at  yield 
(%) 

Stress  at 
break 
(MPa) 

Elongation 
at  break 
(%) 

Young’s 

modulus 

(GPa) 

PYSD 

(%) 

26.4 

NA 

NA 

27.26 

4.6 

1.00 

NA 

(2.06) 

(1.0) 

(0.06) 

43.6 

75.54 

7.2 

51.91 

25.9 

1.71 

32.4 

(0.20) 

(0.1) 

(0.95) 

(12.8) 

(0.09) 

(1.9) 

Values 

in  parentheses  represent  standard  deviation. 

>le  4.13  Results  of  tensile  test  for  polycarbonate  with  chroman-I  end  groups  in  the 

d state 

Mw 

Stress 

Elongation 

Stress  at 

Elongation 

Young’s 

PYSD 

(kg/mole)  at  yield 

at  yield 

break 

at  break 

modulus 

(%) 

(MPa) 

(%) 

(MPa) 

(%) 

(GPa) 

25.9 

NA 

NA 

21.87 

3.4 

1.22 

NA 

(1.30) 

(1.3) 

(0.15) 

44.3 

74.73 

7.5 

50.68 

19.3 

1.65 

32.0 

(0.18) 

(0.0) 

(0.11) 

(4.0) 

(0.05) 

(0.0) 

Values  in  parentheses  represent  standard  deviation. 


greater  density.  After  yielding,  plastic  deformation  cannot  occur  as  much  as  in  the 
quenched  state,  resulting  in  lower  elongation  at  break.  PYSD  in  polycarbonate  results 
from  strain-softening  phenomena  when  samples  start  necking  during  tensile  tests.  Strain- 
softening induces  localized  plastic  strain  instead  of  uniformly  distributed  strain  in  the 
sample.  After  aging,  the  increase  in  PYSD  or  strain  softening  induces  more  localized 
plastic  strain,  which  reduces  the  volume  of  plastic  deformation  in  the  sample.  Hence, 
they  become  more  brittle.  Adam  et  al.  (1975)  proved  that  the  increase  in  PYSD  due  to 
aging  is  not  influenced  by  geometrical  effects  by  measuring  the  draw  ratio  of  the  aged 
and  unaged  samples,  which  come  out  the  same.  Bubeck  et  al.  (1984)  suggested  that  an 
increase  in  PYSD  is  the  result  of  an  increase  in  some  conformations  or  local  molecular 
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packing  arrangements  that  are  unfavorable  for  ductile  deformation.  However,  they  did 
not  demonstrate  proof  of  the  change  in  conformation  they  described.  They  also 
correlated  PYSD  magnitude  with  the  magnitude  of  tan  8 at  80°C.  Lower  PYSD  results  in 
higher  magnitude  of  tan  5 at  80°C  and  higher  resistance  to  aging. 

On  the  other  hand,  samples  with  low  molar  mass,  which  already  exhibit  brittle 
behavior  in  the  quenched  state,  show  no  significant  difference  in  tensile  properties. 

Adam  et  al.  (1975)  compared  the  impact  strength  of  polycarbonate  with  different 
thicknesses.  The  results  show  that  the  samples  in  plane  strain  condition  (sharp  notched, 
high  thickness)  do  not  show  any  difference  in  impact  strength  after  aging  at  130°C  for  48 
hours.  The  samples  in  plane  stress  condition  show  a dramatic  decrease  in  impact 
strength.  However,  the  low  molar  mass  samples  in  this  study  also  become  denser  due  to 
aging  as  shown  in  the  density  results.  At  this  period  of  aging,  the  effect  may  be  too  little 
to  observe.  The  brittle  behavior  is  mainly  caused  by  low  entanglement  density  in  the 
system.  If  the  samples  are  aged  longer,  they  may  show  an  obvious  decrease  in 
elongation. 

Fracture  Toughness  Results 

Fracture  toughness  (KIC)  was  determined  using  quantitative  fractography.  The 
calculation  is  based  on  equations  (2.23)  and  (2.24)  for  polycarbonate  with  low  molar 
mass  and  high  molar  mass,  respectively.  As  expected,  fracture  toughness  of 
polycarbonate  with  high  molar  mass  is  higher  than  polycarbonate  with  low  molar  mass 
for  both  end  groups  (Table  4.14).  To  study  the  effect  of  thermal  history  on  fracture 
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toughness,  the  fracture  toughness  of  polycarbonate  aged  at  Tg-20°C  for  7 days  was  also 
determined  (Table  4.15). 


Table  4.14  Fracture  toughness  of  polycarbonate  in  the  quenched  state 


Sample  code 

End  groups  Mw  (kg/mole) 

KIC  (MPa  ml/2) 

92185 

paracumylphenol 

26.4 

0.78+0.11 

92178 

paracumylphenol 

43.6 

2.83+0.25 

22152 

chroman-I 

25.9 

0.64+0.02 

22141 

chroman-I 

44.3 

3.28+0.32 

1 5 Fracture  toughness  of  polycarbonate 

in  the  aged  state 

Sample  code 

End  groups  Mw  (kg/mole) 

KIC  (MPa  m1/2) 

92185 

paracumylphenol 

26.4 

0.65+0.07 

92178 

paracumylphenol 

43.6 

2.78+0.12 

22152 

chroman-I 

25.9 

0.54+0.04 

22141 

chroman-I 

44.3 

2.54+0.50 

After  aging,  fracture  toughness  tends  to  decrease  for  both  end  groups  in  high 
molar  mass  polycarbonate.  The  change  is  not  so  obvious  in  low  molar  mass 
polycarbonate.  Hill  et  al.  (1990b)  determined  the  K,c  of  commercial  polycarbonate  by 
using  compact  tension  specimen  of  unaged  samples  and  samples  aged  at  120°C  for  240 
hours.  The  values  are  3.23  and  2.63  MPa  ml/2,  respectively.  The  KIC  of  high  molar  mass 
polycarbonate  determined  by  quantitative  fractography  falls  within  the  same  range. 

There  is  no  significant  difference  between  the  two  end  groups  at  the  same  range  of  molar 
mass.  After  cracks  initiate,  they  will  propagate  through  the  sample  until  it  fractures.  For 
brittle  materials,  fracture  occurs  very  fast  due  to  very  little  plastic  deformation  ahead  of 
the  crack  tip.  For  ductile  materials,  plastic  deformation  ahead  of  the  crack  tip  can  occur 
to  the  extent  that  cracks  may  propagate  more  slowly  or  arrest.  Based  on  fracture 
mechanics,  the  size  of  plastic  deformation  ahead  of  the  crack  tip  may  be  calculated.  Irwin 
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suggested  that  the  plastic  deformation  ahead  of  the  crack  tip  can  be  modeled  as  a circle. 
For  plane  strain  condition,  the  radius  of  the  plastic  zone  can  be  calculated  as  follows 
(Kinloch  and  Young,  1983): 


r,  =(1/6^-)(/://c7v)2...(4.3) 

where  r-y  is  the  plastic  zone  radius,  Kx  is  the  stress-intensity  factor,  and  ay  is  yield  stress. 
However,  the  plastic  zone  ahead  of  the  crack  tip  in  many  polymers  is  modeled  as  a line 
rather  than  a circle.  The  calculation  is  based  on  Dugdale’s  analysis  as  follows  (Kinloch 
and  Young,  1983): 

R = (n /%)(Kj  / <jp)2 ...(A A) 

where  R is  the  length  of  plastic  zone,  Kx  is  the  stress-intensity  factor,  and  Op  is  the  internal 
stress,  which  acts  on  the  boundary  of  the  plastic  zone  and  is  assumed  to  be  equal  to  yield 
stress.  The  results  of  the  calculation  for  both  models  are  shown  in  Table  4.16. 


Table  4.16  Size  of  plastic  deformation  at  the  crack  tip  of  polycarbonate  with  yielding  in 
the  quenched  state  and  in  the  aged  state 


Sample  code 

End  groups 

Thermal 

history 

Mw 

(kg/mole) 

ry  (Em) 

R (pm) 

92178 

paracumylphenol 

quenched 

43.6 

112  ±23 

832 ± 173 

92178 

paracumylphenol 

aged 

43.6 

72  ±6 

532  ± 47 

22141 

chroman-I 

quenched 

44.3 

160  ±37 

1182  ± 271 

22141 

chroman-I 

aged 

44.3 

63  ±26 

465  ±192 
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As  shown  in  the  results,  the  plastic  deformation  is  larger  in  the  quenched  state. 
Thermal  history  influences  fracture  behavior,  while  the  influence  of  end  group  is  not 
significant. 

Results  of  Fractal  Dimension  Analysis 

Fracture  surfaces  of  polycarbonate  both  in  quenched  and  aged  states  were 
observed  by  using  AFM.  Fractal  dimensions,  which  are  calculated  by  the  variation 
method  and  the  slit-island  method,  as  well  as  roughness  values  of  the  surfaces,  are  shown 
in  Tables  4.17,  4.18,  4.19,  and  4.20. 


Table  4.17  Fractal  dimension  and  roughness  of  polycarbonate  with  paracumylphenol  end 
groups  in  the  quenched  state 


Mw  Fractal 

(kg/mole)  Dimension* 

Fractal 

Dimensionf 

Roughness1 
RMS  (nm) 

Mean 

Roughness8 

(nm) 

Position  of  the  tip 

26.4 

2.09±0.02 

1.07±0.02 

33±22 

33±22 

in  the  mirror  region 

26.4 

2.16±0.05 

1.10±0.02 

7±3 

7±3 

in  the  mist  region 

43.6 

2.08±0.03 

1.05±0.03 

28±26 

28±26 

0.5  mm  from  the 

originating  flaw 

43.6 

2.17±0.02 

1.10±0.05 

5±0.1 

5±0.1 

Beyond  the  boundary 

that  cracks  change 

plane 

* calculated  based  on  the  variation  method 

+ calculated  based  on  the  slit-island  method 
1 calculated  based  on  equation  3.1 
8 calculated  based  on  equation  3.2 
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Table  4.18  Fractal  dimension  and  roughness  of  polycarbonate  with  chroman-I  end  groups 
in  the  quenched  state 


Mw  Fractal 

(kg/mole)  Dimension* 

Fractal 

Dimension1 

Roughness1 
RMS  (nm) 

Mean 

Roughness5 

(nm) 

Position  of  the  tip 

25.9 

2.14±0.01 

1.14+0.05 

6±1 

6±1 

in  the  mirror  region 

25.9 

2.08+0.01 

1.07±  0.03 

25±7 

25±7 

in  the  mist  region 

44.3 

2.10±0.03 

1.06+  0.03 

22±1 1 

22±1 1 

0.5  mm  from  the 
originating  flaw 

44.3 

* 

2.07±  0.01 

1.02 

14±  4 

19±  5 

Beyond  the 
boundary  that 
cracks  change 
plane 

* calculated  based  on  the  variation  method 


+ calculated  based  on  the  slit-island  method 
1 calculated  based  on  equation  3.1 
5 calculated  based  on  equation  3.2 


Table  4.19  Fractal  dimension  and  roughness  of  polycarbonate  with  paracumylphenol  end 
groups  in  the  aged  state 


Mw  Fractal  Fractal  Roughness1  Mean 

(kg/mole)  Dimension*  Dimension1  RMS  (nm)  Roughness* * * 5 


(nm) 


26.4 

2.12+0.05 

1.07+  0.04 

17+15 

20+18 

26.4 

2.15+0.00 

1.09+0.02 

13+7 

17+10 

43.6 

2.20+0.07 

1.17+0.05 

13+1 

14+4 

43.6 

2.12+0.01 

1.14+0.03 

11+3 

14+4 

* calculated  based  on  the  variation  method 

+ calculated  based  on  the  slit-island  method 

1 calculated  based  on  equation  3.1 

5 calculated  based  on  equation  3.2 


Position  of  the  tip 


in  the  mirror  region 
in  the  mist  region 
0.5  mm  from  the 
originating  flaw 
Beyond  the 
boundary  that 
cracks  change 
plane 
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Table  4.20  Fractal  dimension  and  roughness  of  polycarbonate  with  chroman-I  end  groups 
in  the  aged  state 


Mw  Fractal 

(kg/mole)  Dimension* 

Fractal 

Dimension1' 

Roughness* 
RMS  (nm) 

Mean 

Roughness8 

(nm) 

Position  of  the  tip 

25.9 

2.10±0.02 

1.07±0.02 

3±1 

4±1 

in  the  mirror  region 

25.9 

2.08±0.03 

1.10±0.03 

58±15 

74±23 

in  the  mist  region 

44.3 

2.12±0.02 

1.08±0.01 

1±0.3 

2±0.3 

0.5  mm  from  the 
originating  flaw 

44.3 

2.18±0.03 

1.16±0.04 

3±0.3 

4±0.4 

Beyond  the 
boundary  that  cracks 
change  plane 

* calculated  based  on  the  variation  method 


+ calculated  based  on  the  slit-island  method 
* calculated  based  on  equation  3.1 
5 calculated  based  on  equation  3.2 


Both  roughness  values  depend  on  the  scale  of  the  measurement.  Different 
sampling  lengths  may  yield  different  values  for  the  same  profile  or  surface.  Figure  4.20 
shows  a single  simulated  Brownian  motion  profile  with  two  different  RMS  values 
depending  on  the  scale  of  the  measurement.  However,  fractal  dimension  remains 


constant.  Another  limitation  of  RMS  is  shown  in  Figure  4.21 . Three  different  profiles 
were  created,  each  having  three  different  frequency  components  (k)  and  amplitudes  (A), 
in  which  A,=T00,  A2=50,  and  A3=25.  They  all  have  the  same  RMS  value  although  they 
are  totally  different.  Clearly,  the  RMS  parameter  fails  to  distinguish  certain  spatial 
frequency  differences.  As  described  above,  fractal  dimension  represents  the  average 
property  of  the  entire  fracture  surface  and  measures  its  tortuosity  better  than  surface 
roughness  does. 


The  results  in  each  category  are  obtained  from  two  samples  and  several  scanning 
areas  in  each  sample.  There  is  a large  variation  in  roughness  values  even  at  the  same 
distance  from  the  originating  flaw,  depending  on  the  position  of  the  tip.  Fractal 
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Figure  4.20  Roughness  value  and  fractal  dimension  of  a single  simulated  Brownian 
motion  profile.  Source:  Spanos  and  Irene  (1993) 
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Figure  4.21  RMS  values  of  three  generated  profiles  with  different  frequency  components 
and  amplitudes.  Source:  Spanos  and  Irene  (1993) 
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dimension  plots  for  all  fractal  dimension  values  are  straight  lines,  as  in  the  examples 
shown  in  Figures  3.4  and  3.5.  Fractal  dimension  from  the  variation  method  and  the  slit- 
island  method  can  be  compared  by  looking  at  the  fractional  part  of  the  fractal  dimension, 
e.g.,  fractal  dimension  2.13  from  the  variation  method  is  equal  to  fractal  dimension  1.13 
from  the  slit-island  method.  Most  of  the  fractal  dimensions  from  the  slit-island  method  in 
this  study  are  lower  than  those  from  the  variation  method  but  they  follow  the  same 
patterns.  The  calculation  by  the  variation  method  is  based  on  the  volume  occupied  by  the 
real  surface,  while  the  slit-island  is  based  on  the  fractal  dimension  of  the  contour  line 
after  slicing  the  fracture  surface. 

Fractal  dimension  values  vary  with  distance  from  the  originating  flaw,  which  may 
be  due  to  the  difference  in  crack  velocity  as  it  propagates  through  the  sample  or  due  to 
data  scattering.  However,  the  propagation  of  the  crack  should  contribute  to  the  tortuosity 
on  the  fracture  surface. 

The  Relationship  between  Fractal  Dimension  and  Fracture  Toughness 

As  described  in  the  background  chapter,  there  is  still  some  controversy  regarding 
the  relationship  between  fractal  dimension  and  fracture  toughness  for  polymers.  By  using 
the  average  values  of  the  fractal  dimension  from  different  regions  in  each  type  of  sample, 
the  plots  between  (D*)1'2  (the  fractional  part  of  the  fractal  dimension)  and  KIC  for  both 
types  of  fracture  can  be  determined.  These  plots  are  shown  in  Figures  4.22  and  4.23. 

After  aging,  fractal  dimensions  of  polycarbonate  with  high  molar  mass  tend  to 
increase.  On  the  other  hand,  fractal  dimensions  of  polycarbonate  with  low  molar  mass 
tend  to  decrease  after  aging.  There  are  not  enough  data  points  to  fit  the  equation  or  to  do 
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(D*)1/2 


Figure  4.22  Plot  between  (D*)l/2  and  KIC  of  polycarbonate  with  paracumylphenol  end 
groups,  Mw=43.6  kg/mole  and  polycarbonate  with  chroman-I  end  groups,  Mw=44.3 
kg/mole  in  the  quenched  and  aged  states. 
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Figure  4.23  Plot  between  (D*)1/2  and  KIC  of  polycarbonate  with  paracumylphenol  end 
groups,  Mw=26.4  kg/mole  and  polycarbonate  with  chroman-I  end  groups,  Mw=25.9 
kg/mole  in  the  quenched  and  aged  states. 
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the  regression  analysis.  The  relationship  between  fractal  dimension  and  fracture 
toughness  of  polycarbonate  with  high  molar  mass,  which  have  ductile  fracture,  exhibits  a 
negative  trend,  as  found  in  some  metals.  The  relationship  between  fractal  dimension  and 
fracture  toughness  of  polycarbonate  with  low  molar  mass,  which  have  brittle  fracture, 
demonstrates  a positive  trend,  as  found  in  most  glass  and  glass-ceramics.  Therefore, 
fractal  dimension  is  based  on  the  fracture  mechanism  more  than  on  materials. 

Fractography  Results 

Fractography  at  the  Macroscopic  Level 

The  optical  micrographs  of  fracture  surfaces  of  polycarbonate  with  both  end 
groups  in  the  quenched  and  aged  states  are  shown  in  Figures  4.24,  4.25,  4.26,  4.27,  4.28, 
4.29,  4.30,  and  4.3 1 . The  micrographs  show  two  different  patterns  of  fracture,  brittle  and 
ductile.  Brittle  fracture  surfaces  consist  of  the  originating  flaw,  the  mirror  or  smooth  area, 
the  mist  or  the  area  with  fine  markings,  and  the  hackle  or  rough  area.  Polycarbonate  with 
low  molar  mass  exhibits  this  type  of  fracture  both  in  the  quenched  and  aged  state.  There 
is  no  difference  shown  between  end  groups  in  the  fracture  pattern.  Ductile  fracture 
surfaces  consist  of  the  originating  flaw,  the  traces  of  the  cracks  that  propagate  from  the 
originating  flaw  and  grow  larger,  and  the  abrupt  transition  before  the  cracks  change  the 
plane.  Polycarbonate  with  high  molar  mass  exhibits  this  type  of  fracture  both  in  the 
quenched  and  aged  states. 

Fracture  surfaces  of  polycarbonates  with  brittle  fracture  exhibit  features  similar  to 
those  in  inorganic  glass.  Cracks  initiate  from  pre-existing  flaws  and  propagate  through 
the  mirror,  mist,  and  hackle  regions  until  catastrophic  fracture  occurs.  The  velocity  of 
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Figure  4.24  Fracture  surface  of  polycarbonate  with  paracumylphenol  end  groups 
Mw=26.4  kg/mole  in  the  quenched  state.  (Optical  micrograph  at  50x) 
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Figure  4.25  Fracture  surface  of  polycarbonate  with  paracumylphenol  end  groups 
Mw=43.6  kg/mole  in  the  quenched  state.  (Optical  micrograph  at  50x) 
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Figure  4.26  Fracture  surface  of  polycarbonate  with  chroman-I  end  groups  Mw=25.9 
kg/mole  in  the  quenched  state.  (Optical  micrograph  at  50x) 
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Figure  4.27  Fracture  surface  of  polycarbonate  with  chroman-I  end  groups  Mw=44.3 
kg/mole  in  the  quenched  state.  (Optical  micrograph  at  50x) 
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Figure  4.28  Fracture  surface  of  polycarbonate  with  paracumylphenol  end  groups 
Mw=26.4  kg/mole  in  the  aged  state.  (Optical  micrograph  at  50x) 
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Figure  4.29  Fracture  surface  of  polycarbonate  with  paracumylphenol  end  groups 
Mw=43.6  kg/mole  in  the  aged  state.  (Optical  micrograph  at  50x) 
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Figure  4.30  Fracture  surface  of  polycarbonate  with  chroman-I  end  groups  Mw=25.9 
kg/mole  in  the  aged  state.  (Optical  micrograph  at  50x) 
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Figure  4.31  Fracture  surface  of  polycarbonate  with  chroman-I  end  groups  Mw=44.3 
kg/mole  in  the  aged  state.  (Optical  micrograph  at  50x) 
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the  crack  increases  from  the  mirror  to  the  mist  to  the  hackle  regions.  In  the  hackle  region, 
cracks  usually  branch  and  propagate  in  different  planes  which  can  be  seen  from  the  rough 
surface.  The  size  of  the  mirror  region  is  small  compared  to  the  mist  or  hackle  region. 

Fracture  surfaces  of  polycarbonate  with  ductile  fracture  show  distinctive  features 
from  brittle  fracture.  The  originating  flaw,  which  may  be  inherent  or  may  be  caused  by 
processing,  is  usually  similar  in  size  to  that  in  brittle  fracture.  However,  cracks  propagate 
through  the  area  surrounding  the  flaw  to  a larger  size  with  relative  low  velocity,  then 
change  plane  and  increase  in  velocity  after  the  boundary.  During  necking,  several  cracks 
emanate  from  the  slow  crack  growth  region  until  they  reach  the  boundary,  and  then  they 
change  to  fast  crack  growth  until  fracture  (Sherman  et  al.,  1982).  Figure  4.32  shows  the 
surface  beyond  the  boundary  that  cracks  change  plane  and  velocity.  The  change  of  plane 
of  propagation  and  crack  branching  can  be  seen  clearly,  which  indicates  that  the  crack 
velocity  is  high  compared  to  the  region  before  the  boundary. 

From  the  observation  of  fracture  surfaces  at  the  macroscopic  level,  the  aging 
condition  in  this  study  does  not  alter  the  fracture  mechanism  of  the  samples.  Low  molar 
mass  polycarbonate  exhibits  brittle  fracture  both  in  quenched  and  aged  states.  High 
molar  mass  polycarbonate  exhibits  ductile  fracture  both  in  quenched  and  aged  state; 
however,  they  become  more  brittle  as  shown  in  tensile  test  results  and  fracture  toughness 
results.  The  longer  aging  time  may  make  the  samples  more  brittle  until  the  brittle 


fracture  can  be  observed  from  the  fracture  surfaces. 


128 


Figure  4.32  Fracture  surface  of  polycarbonate  with  chroman-I  end  groups  Mw=44.3 
kg/mole  in  the  aged  state  beyond  the  boundary  that  cracks  change  plane.  (Optical 
micrograph  at  5 Ox) 
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Fractography  at  the  Molecular  Level 

AFM  images  of  polycarbonates  with  brittle  and  ductile  fracture  in  the  quenched 
state  are  shown  in  Figures  4.33  — 4.40  and  AFM  images  of  polycarbonates  with  brittle 
and  ductile  fracture  in  the  aged  state  are  shown  in  Figures  4.41  - 4.48.  The  images 
selected  are  representative  of  the  images  scanned  at  the  same  distance  from  the 
originating  flaw.  These  images  confirm  the  irregular  nature  of  the  fracture  surfaces  even 
at  the  nanometer  scale  range.  The  hypothesis  that  the  fracture  process  has  fractal 
characteristics  is  reasonable.  The  position  of  the  scanning  area  can  make  the  images  look 
very  different,  depending  on  the  surface  features  that  the  tip  scanned.  Characterization  of 
fracture  surfaces  using  only  the  images  may  be  difficult.  Therefore,  quantitative 
parameters  such  as  fractal  dimension  or  roughness  are  rather  used.  However,  the  images 
alone  can  be  a useful  guideline  in  comparing  different  materials  or  the  same  materials 
under  different  conditions. 

Fractography  of  brittle  polycarbonate  exhibits  the  features  similar  to  those  of 
soda-lime  glass  studied  by  Kulawansa  et  al.  (1994).  The  ridges  and  steps  are  the 
prominent  features.  They  suggested  that  crack  fingering  is  responsible  for  these  features. 
Crack  fingering  results  from  crack  front  instability,  which  occurs  when  there  is  a positive 
stress  gradient.  The  crack  develops  the  fingers,  which  move  towards  the  high  stress 
region.  The  ridges  represent  these  crack  fingers.  Changes  in  the  direction  of  maximum 
principal  stress  rotate  the  planes  of  the  crack  fingers.  After  the  maximum  principal  stress 
is  returned  to  the  normal  direction,  these  crack  fingers  are  not  coplanar  any  more, 
producing  step  features.  The  diagram  of  step  formation  is  shown  in  Figure  4.49. 
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Figure  4.33  AFM  image  of  fracture  surface  of  polycarbonate  with  paracumylphenol  end 
groups  Mw=26.4  kg/mole  in  the  quenched  state,  scanning  position  is  in  the  mirror  region. 
The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.34  AFM  image  of  fracture  surface  of  polycarbonate  with  paracumylphenol  end 
groups  Mw=26.4  kg/mole  in  the  quenched  state,  scanning  position  is  in  the  mist  region. 
The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.35  AFM  image  of  fracture  surface  of  polycarbonate  with  paracumylphenol  end 
groups  Mw=43.6  kg/mole  in  the  quenched  state,  scanning  position  ~ 0.5  mm  from  the 
originating  flaw.  The  arrow  indicates  the  general  direction  of  crack  growth. 


133 


nM 


X 200.000  nM/diu 
Z 50.000  nM/diu 


92178q2*.4 


Figure  4.36  AFM  image  of  fracture  surface  of  polycarbonate  with  paracumylphenol  end 
groups  Mw=43.6  kg/mole  in  the  quenched  state,  scanning  position  is  beyond  the  boundary 
that  cracks  change  plane.  The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.37  AFM  image  of  fracture  surface  of  polycarbonate  with  chroman-I  end  groups 
Mw=25.9  kg/mole  in  the  quenched  state,  scanning  position  is  in  the  mirror  region.  The 
arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.38  AFM  image  of  fracture  surface  of  polycarbonate  with  chroman-I  end  groups 
Mw=25.9  kg/mole  in  the  quenched  state,  scanning  position  is  in  the  mist  region.  The 
arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.39  AFM  image  of  fracture  surface  of  polycarbonate  with  chroman-I  end  groups 
Mw=44.3  kg/mole  in  the  quenched  state,  scanning  position  ~ 0.5  mm  from  the  originating 
flaw.  The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.40  AFM  image  of  fracture  surface  of  polycarbonate  with  chroman-I  end  groups 
Mw=44.3  kg/mole  in  the  quenched  state,  scanning  position  is  beyond  the  boundary  that 
cracks  change  plane.  The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.41  AFM  image  of  fracture  surface  of  polycarbonate  with  paracumylphenol  end 
groups  Mw=26.4  kg/mole  in  the  aged  state,  scanning  position  is  in  the  mirror  region.  The 
arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.42  AFM  image  of  fracture  surface  of  polycarbonate  with  paracumylphenol  end 
groups  Mw=26.4  kg/mole  in  the  aged  state,  scanning  position  is  in  the  mist  region.  The 
arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.43  AFM  image  of  fracture  surface  of  polycarbonate  with  paracumylphenol  end 
groups  Mw=43.6  kg/mole  in  the  aged  state,  scanning  position  ~ 0.5  mm  from  the 
originating  flaw.  The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.44  AFM  image  of  fracture  surface  of  polycarbonate  with  paracumylphenol  end 
groups  Mw=43.6  kg/mole  in  the  aged  state,  scanning  position  is  beyond  the  boundary  that 
cracks  change  plane.  The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.45  AFM  image  of  fracture  surface  of  polycarbonate  with  chroman-I  end  groups 
Mw=25.9  kg/mole  in  the  aged  state,  scanning  position  is  in  the  mirror  region.  The  arrow 
indicates  the  general  direction  of  crack  growth. 
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2 200.000  nM/diu 
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Figure  4.46  AFM  image  of  fracture  surface  of  polycarbonate  with  chroman-I  end  groups 
Mw=25.9  kg/mole  in  the  aged  state,  scanning  position  is  in  the  mist  region.  The  arrow 
indicates  the  general  direction  of  crack  growth. 
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Figure  4.47  AFM  image  of  fracture  surface  of  polycarbonate  with  chroman-I  end  groups 
Mw=44.3  kg/mole  in  the  aged  state,  scanning  position  ~0.5  mm  from  the  originating  flaw. 
The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.48  AFM  image  of  fracture  surface  of  polycarbonate  with  chroman-I  end  groups 
Mw=44.3  kg/mole  in  the  aged  state,  scanning  position  is  beyond  the  boundary  that  cracks 
change  plane.  The  arrow  indicates  the  general  direction  of  crack  growth. 
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Figure  4.49  Diagram  of  step  formation.  Source:  Kulawansa  et  al.  (1994) 
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Fractography  of  ductile  polycarbonate  displays  the  similar  features  of  ridges  and 
steps  as  in  brittle  polycarbonate,  but  it  usually  includes  the  bump  features.  The 
relationship  between  bumps  and  the  fracture  process  is  unclear.  Kulawansa  et  al.  (1992) 
suggested  that  they  also  result  from  crack  instability. 

The  difference  between  the  features  in  the  quenched  and  the  aged  states  of  ductile 
polycarbonate  is  interesting.  Figures  4.35  and  4.39,  which  are  polycarbonate  in  the 
quenched  state,  have  less  tortuosity  than  those  in  Figures  4.43  and  4.47,  which  are  in  the 
aged  state.  Figures  4.43  and  4.47  show  tortuosity  on  a very  small  scale.  The  higher  yield 
stress  in  the  aged  state  may  contribute  to  these  tortuosities.  The  chains  stretched  more 
before  the  polymer  started  necking,  resulting  in  more  numerous  deformed  marks  on  the 
surfaces. 


CHAPTER  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  completion  of  this  dissertation  has  accomplished  two  goals.  The  first  goal 
was  to  study  ductile-to-brittle  transition  in  amorphous  thermoplastics  due  to  internal 
variables,  i.e.,  molar  mass,  and  external  variables,  i.e.,  thermal  history.  A ductile-to- 
brittle  transition  due  to  molar  mass  is  prominent  in  tensile  properties,  which  change  from 
failure  with  yielding,  necking,  and  high  elongation  to  failure  without  yielding,  necking, 
and  low  elongation.  Fracture  toughness  also  decreases  as  a function  of  molar  mass. 
Fractography  at  the  macroscopic  level  changes  from  the  pattern  of  ductile  fracture,  in 
which  cracks  propagate  to  a larger  size  before  they  change  direction  and  fracture  occurs, 
to  the  pattern  of  brittle  fracture  with  mirror,  mist,  hackle  regions.  Fractography  at  the 
molecular  level  of  polycarbonate  with  high  molar  mass  has  less  tortuosity  than 
polycarbonate  with  low  molar  mass  at  the  same  distance  from  the  flaw.  This  may  be  due 
to  lower  crack  velocity  in  polycarbonate  with  ductile  fracture  than  in  polycarbonate  with 
brittle  fracture. 

A ductile-to-brittle  transition  due  to  thermal  history  or  physical  aging  in  this  study 
is  prominent  only  in  polycarbonate  with  high  molar  mass.  Increase  in  yield  stress, 
decrease  in  elongation,  and  fracture  toughness  are  the  main  changes.  Polycarbonate  with 
low  molar  mass  does  not  exhibit  significant  changes  in  tensile  properties  or  fracture 
toughness,  but  there  is  a trend  of  decreasing  elongation  and  fracture  toughness. 
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Fractography  at  the  macroscopic  level  does  not  exhibit  significant  differences 
between  polycarbonate  in  the  quenched  state  and  the  aged  state.  Fractography  of  low 
molar  mass  polycarbonate  exhibits  a brittle  pattern  with  mirror,  mist,  hackle  regions  both 
in  the  quenched  state  and  in  the  aged  state.  Fractography  of  high  molar  mass 
polycarbonate  exhibit  a ductile  pattern,  in  which  cracks  propagate  to  a larger  size  before 
they  change  direction  and  fracture  occurs,  both  in  the  quenched  state  and  in  the  aged 
state.  However,  fractography  at  the  molecular  level  of  polycabonate  with  high  molar 
mass  exhibits  higher  tortuosity  in  the  aged  state  than  in  the  quenched  state.  This  may  be 
due  to  higher  yield  stress  in  the  aged  state,  which  makes  the  chains  stretch  more  before 
the  polymer  starts  necking,  resulting  in  more  numerous  deformed  marks  on  the  surfaces. 

Fractography  at  the  molecular  level  reveals  that  the  fracture  process  has  fractal 
characteristics.  The  irregular  topography  at  the  nanometer  scale  can  be  seen  in  the  area 
that  looks  smooth  macroscopically.  Crack  propagation  induces  the  features  or  markings 
on  the  fracture  surfaces  even  at  the  molecular  level.  Atomic  force  microscopy  is  a 
valuable  tool  for  polymers  since  surface  coating  is  not  required.  Tappingmode®  AFM, 
which  applies  only  nanoNewton  force  on  the  surface,  is  highly  appropriate  for  studying 
fracture  surfaces.  The  tip  does  not  damage  the  features  of  the  fracture  surfaces. 

However,  since  fracture  surfaces  are  usually  not  in  the  same  plane,  the  positioning  of  the 
tip  is  quite  difficult. 

Fractal  dimension  is  a quantitative  parameter  used  to  characterize  fracture 
surfaces.  It  is  best  to  have  a quantitative  parameter  paired  with  visual  observation.  The 
relationship  between  fractal  dimension  and  fracture  toughness  has  a positive  trend  in 
polycarbonate  with  brittle  fracture  but  a negative  trend  in  polycarbonate  with  ductile 
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fracture.  Fractal  dimensions  relate  to  fracture  toughness  differently,  depending  on  the 
fracture  mechanism.  However,  the  data  points  are  too  few  to  fit  the  equation  or  to  do  the 
regression  analysis.  The  experiment  with  more  samples  with  different  fracture  toughness 
values  should  be  done  in  order  to  develop  the  mathematical  model  for  the  relationship. 

Quantitative  fractography,  which  uses  the  relationship  between  fracture  stress  and 
flaw  size,  provides  reliable  fracture  toughness  values  both  for  brittle  and  ductile 
polycarbonates.  This  technique  has  not  been  reported  with  ductile  polymers  before.  The 
results  are  comparable  to  the  fracture  toughness  values  in  the  literature  determined  by 
inserting  the  sharp  crack  in  the  sample.  Quantitative  fractography  can  be  done  more 
easily  since  no  special  geometry  of  the  samples  is  needed. 

The  second  goal  that  this  dissertation  has  accomplished  is  to  study  the  influence 
of  end  groups  on  physical  properties,  mechanical  properties,  and  the  aging  rate.  There  is 
the  important  influence  of  end  groups  on  the  temperature  and  the  activation  energy  at  the 
glass  transition,  density,  and  the  shape  and  the  activation  energy  of  P-relaxation.  Large 
end  groups  induce  high  glass  transition  with  high  activation  energy  but  broader  p- 
relaxation  with  low  activation  energy.  Large  end  groups  also  induce  lower  density  upon 
quenching  from  above  the  glass  transition.  However,  tensile  properties,  fracture 
toughness,  fractal  dimension,  and  fractography  at  the  macroscopic  level  and  the 
molecular  level  are  not  significantly  different  between  polycarbonate  with  two  end 
groups  used  in  the  study  at  the  same  range  of  molar  mass.  The  fact  that  end  groups  have 
no  influence  on  properties  described  above  is  desirable  for  using  the  end  group  as  a factor 
to  control  physical  aging  process  by  retaining  the  other  major  properties.  Unfortunately, 
the  measured  rates  of  physical  aging  of  polycarbonate  with  two  different  end  groups  used 
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in  this  study  do  not  show  any  significant  difference.  The  sizes  of  the  two  end  groups 
used  in  this  study  may  not  be  different  enough  from  each  other  to  cause  a difference  in 
the  rate  of  physical  aging.  Experiments  with  larger  end  groups  need  to  be  done  to 
determine  their  effect  on  the  rate  of  physical  aging. 


APPENDIX  A 

RAW  DATA  OF  FRACTURE  TOUGHNESS  DETERMINATION 


Fracture  surface  analysis  by  optical  microscopy  of  polycarbonates  with  brittle  fracture  in  the  quenched  and  aged  state 
Sample  End  Flaw  size  R mirror  Rmist  R hackle  Flaw  type  KIC(brittle) 
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Remarks: 


sample 

q = quenched  state 

a = aged  state 

end  group 

P = paracumylphenol 

C = Chroman-I 

Type  of  flaw 

i = inclusion 

c = crack 

b = bubble 

r'y  = Irwin’s  plastic  zone  radius 

R = Dugdale’s  plastic  zone  length 


APPENDIX  B 

RAW  DATA  OF  FRACTAL  DIMENSION  AND  ROUGHNESS  VALUES 


Fractal  dimension  and  roughness  of  polycarbonates  with  brittle  fracture  in  the  quenched  and  aged  state 

Image  file  Fractal  Average  Fractal  Average  Mean  RMS4  Position  of  scanning 

name  dimension1  fractal  dimension2  fractal  roughness3  (RJ  (nm) 
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Image  file  Fractal  Average  Fractal  Average  Mean  RMS4  Position  of  scanning 

name  dimension1  fractal  dimension2  fractal  roughness3  (Rq)  (nm) 


158 


a 

c 

a 

o 

o 

o 

c 

c 

c 

c 

'5b 

'3) 

‘5b 

o 

o 

o 

• ^4 

.2 

<U 

u 

<D 

00 

*5b 

CJ) 

’5b 

S-i 

i-i 

<D 

0) 

<D 

<D 

k4 

t-H 

i-> 

i-i 

Fh 

o 

o 

o 

4-* 

4-* 

4-» 

4-» 

fc 

fc 

fc 

CO 

• 

C/5 

r. 

CO 

£ 

o 

£ 

u 

£ 

£ 

<L> 

£ 

<u 

£ 

tu 

£ 

<L> 

O 

JO 

JO 

JO 

^C3 

JO 

JO 

JO 

4-» 

4-» 

4-> 

4— » 

■*-> 

4— » 

a 

.s 

.s 

S3 

.s 

.s 

r-  m - $ 

OO  O ON  “ 

h in  n 


m m <N 


rn 

NO 

m 


oo 


NO  t"- 
O O 
<N  ON 

K O 
on  in 


CN  a *n  2)  a 


ON 


<N  NO  00 
ij-  cn 


CN 

o 
o 
+ 1 
r- 
o 


m (N  NO 

no  on  m 

o o o 


oo 

NO 

in 

■*t 


in 

o 

NO 


CN 
<N  it 

co 

r~~  tj- 


m 

o 

o 


NO  it 
NO  CN 
O i-i 


m 

o 

o 


+ 1 


oo 

p 

CN 


cn  m it  no  on  —i  r>- 
— - t"  o no  cn  oo  m 
p <-<  o — < o o 
on  cn  (N  (-sj  H on  cn 


>— 1 m 

^ ^ (N  N m rn 

<n  £ it  cn  £ Tt 

B o ts  n nj 

(N  B fS  tS  (N  CB  (N 

in  cn  >n  m in  on  in 

— < in  ^ in  i 

CN  ^ (N  ON  (N  -H  CN 

M tN  M (N  N M M 

CN  CN 


159 


Remarks: 


Image  file  name 

92185  = polycarbonate  with  paracumylphenol  end  group  Mw  26.4  kg/mole 
22152  = polycarbonate  with  chroman-I  end  group  Mw  25.9  kg/mole 
q = quenched  state 
a = aged  state 

Fractal  dimension1  = fractal  dimension  from  the  variation  method 
Fractal  dimension2  = fractal  dimension  from  the  slit-island  method 
Mean  roughness3  = calculated  based  on  the  equation  (3.1) 

RMS4=  calculated  based  on  the  equation  (3.2) 

* = cannot  calculate  fractal  dimension 


Fractal  dimension  and  roughness  of  polycarbonates  with  ductile  fracture  in  the  quenched  and  aged  state 
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Image  file  Fractal  Average  Fractal  Average  Mean  RMS4  Position  of  scanning 

name  dimension'  fractal  dimension2  fractal  roughness3  (Rq)  (nm) 
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Remarks: 


Image  file  name 

92178  = polycarbonate  with  paracumylphenol  end  group  Mw  43.6 
kg/mole 

22141  = polycarbonate  with  chroman-I  end  group  Mw  44.3 
kg/mole 

q = quenched  state 
a = aged  state 

Fractal  dimension'  = fractal  dimension  from  the  variation  method 
Fractal  dimension2  = fractal  dimension  from  the  slit-island  method 
Mean  roughness3  = calculated  based  on  the  equation  (3.1) 

RMS4=  calculated  based  on  the  equation  (3.2) 

* = cannot  calculate  fractal  dimension 


APPENDIX  C 

CHARACTERISTIC  LENGTH  DETERMINATION 
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According  to  the  empirical  relationship  between  fracture  toughness  and  fractal 
dimension  in  the  glass  and  glass  ceramics  developed  by  Mecholsky  et  al.  (1989),  the 
characteristic  length  can  be  determined  from  the  following  equation: 

KIC  =Ea0U2D*'12 

where  KIcis  fracture  toughness  in  MPa  m1/2,  E is  Young’s  modulus  in  MPa,  a0  is  the 
characteristic  length  in  meter,  and  D*  is  the  fractional  part  of  fractal  dimension 

Polycarbonate  with  brittle  fracture  in  the  study  exhibits  similar  relationship 
between  fracture  toughness  and  fractal  dimension  as  glass  and  glass  ceramics.  Hence,  the 
characteristic  length  can  be  calculated  as  shown  in  Table  I. 


Table  I Calculated  a0  for  polycarbonate  with  brittle  fracture 


End 

group 

Molar 

mass 

(kg/mole) 

Thermal 

history 

KIC 

(MPa  ml/2) 

E 

(MPa) 

D* 

a0 

(A) 

p-cumyl 

26.4 

quenched 

0.78 

1 .08*  1 03 

0.13 

40100 

p-cumyl 

26.4 

aged 

0.65 

1 .00*  1 03 

0.13 

32500 

Chroman-I 

25.9 

quenched 

0.64 

1 .06*  1 03 

0.11 

33100 

Chroman-I 

25.9 

aged 

0.54 

1 .22*  103 

0.09 

21800 

The  SEM  micrographs  of  fracture  surfaces  are  shown  in  Figures  1 and  2.  The 
distances  shown  in  Figures  1 and  2,  which  are  in  the  range  of  1000  nm  or  10000  A,  may 
correspond  to  a0.  Fracture  will  take  part  along  the  space  or  along  the  path  that  requires 
the  minimum  energy  in  the  system,  which  can  be  along  the  free  volume  in  amorphous 
thermoplastics.  Another  observation  from  the  calculated  characteristic  length 
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Figure  1 Fracture  surface  close  to  the  initiating  crack  region  of  polycarbonate  with 
paracumylphenol  end  groups,  Mw=26.4  kg/mole  in  the  aged  state.  (SEM  micrograph  at 
30,000x) 
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Figure  2 Fracture  surface  close  to  the  initiating  crack  region  of  polycarbonate  with 
chroman-I  end  groups,  Mw=25.9  kg/mole  in  the  aged  state.  (SEM  micrograph  at  30,000x) 
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is  the  characteristic  length  decreases  after  aging  in  both  end  groups,  which  may  result 
from  a decrease  in  free  volume. 

Polycarbonate  with  ductile  fracture  exhibits  an  inverse  relationship  between 
fracture  toughness  and  fractal  dimension.  Therefore,  the  above  mentioned  equation  is  not 
valid. 


APPENDIX  D 
T-TEST  ANALYSIS 
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In  order  to  determine  whether  there  is  a significant  difference  in  the  results,  t-test 
was  performed  for  density,  fracture  toughness,  fractal  dimension  and  activation  energy 
results.  P values  of  all  t-tests  are  shown  in  Tables  I and  II. 


Table  I P value  for  the  t-test  of  density,  fracture  toughness  and  fractal  dimension  results 


Sample 

Comparison 

P value 
for  density 

P value 
for  KIC 

P value 

for  fractal  dimension 

1 85q/l  85a 

Thermal  history 

0.0135 

0.10 

0.67 

1 52q/l  52a 

9.2619e-5 

0.02 

0.30 

141q/141a 

1.0083e-6 

0.11 

0.03 

1 78q/l  78a 

0.3683 

0.74 

0.18 

185q/152q 

End  groups 

8.6894e-6 

0.14 

0.44 

1 78q/l  4 1 q 

5.5120e-7 

0.19 

0.63 

185a/152a 

1.4943e-3 

0.04 

0.03 

178a/141a 

0.0607 

0.31 

0.62 

185q/178q 

Molar  mass 

2.2143e-3 

4.2367e-4 

0.44 

152q/141q 

1.2479e-3 

3.0822e-4 

0.53 

1 85a/l  78a 

0.5468 

5.0000e-10 

0.67 

152a/141a 

4.9844e-4 

2.3510e-3 

0.02 

185  = polycarbonate  with  paracumylphenol  end  groups,  Mw=  26.4  kg/mole 
178  = polycarbonate  with  paracumylphenol  end  groups,  Mw=  43.6  kg/mole 
152  = polycarbonate  with  chroman-I  end  groups,  Mw=  25.9  kg/mole 
141  = polycarbonate  with  chroman-I  end  groups,  Mw=  44.3  kg/mole 
q = quenched 
a = aged 


Table  II  P value  for  the  t-test  of  activation  energy  results 


Sample 

Comparison 

P value 

183q/128q 

Activation  energy  of  the  a-  relaxation 
as  a function  of  end  groups 

0.05 

183q/128q 

Activation  energy  of  the  P-  relaxation 
as  a function  of  end  groups 

0.03 

1 83  = polycarbonate  with  paracumylphenol  end  groups,  Mw=  3 1 .2  kg/mole 
128  = polycarbonate  with  chroman-I  end  groups,  Mw=  33.2  kg/mole 
q = quenched 
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Based  on  the  value  of  a = 0.05,  there  is  a significant  difference  in  the  density  as  a 
function  of  thermal  history,  end  groups  and  molar  mass  in  all  categories  of  samples 
except  in  high  molar  mass  polycarbonate  with  paracumylphenol  end  groups.  This 
category  exhibits  unusual  behavior  after  aging.  More  samples  may  be  used  to  determine 
whether  their  densities  are  the  same  before  and  after  aging.  There  is  a significant 
difference  in  fracture  toughness  as  a function  of  molar  mass  but  not  as  a function  of 
thermal  history  or  end  groups.  Thermal  history  or  end  groups  used  in  the  study  may  not 
create  significant  changes  in  fracture  toughness.  Modifying  the  aging  condition  or  the 
size  of  end  groups  may  help  verifying  the  discrepancies.  There  is  no  significant 
difference  in  the  fractal  dimension  for  all  categories  of  samples  except  1 4 1 q/ 1 4 1 a, 
185a/152a,  and  152a/141a.  More  samples  and  more  measurements  in  one  sample  should 
be  used.  The  change  in  fractal  dimension  may  be  too  small  for  the  range  of  fracture 
toughness  in  the  study.  Samples  with  wider  range  of  fracture  toughness  should  be  used 
as  well.  There  is  a significant  difference  in  the  activation  energy  of  a and  P relaxation  as 
a function  of  end  groups  (Table  II). 
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